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SUMMARY 


Under  a  tri-service  cooperative  effort  sponsored  by  the  Air 
Force,  Army,  and  Navy  with  extensive  participation  by  NASA,  a 
missile-fin  data  base  for  a  wide  range  of  configurations  and  flow 
conditions  was  obtained.  These  data  were  incorporated  into  an 
engineering  method  for  predicting  the  aerodynamic  characteristics 
of  typical  cruciform  missile  configurations  over  a  wide  range  of 
angles  of  attack,  fin  deflection  angles,  roll  angles,  and  Mach 
numbers.  This  final  report  documents  the  test  programs, 
describes  the  new  code  MISSILE  3,  and  presents  comparisons  of 
independent  experiment  and  predicted  results  to  verify  the 
code.  A  user's  manual  for  the  code  is  included. 


FOREWORD 

This  is  the  final  technical  report  on  the  work  performed 
under  Contract  N00014-80-C-0700  from  16  June  1980  to  30  July 
1986.  Dr.  Robert  Whitehead  of  the  Office  of  Naval  Research  acted 
as  the  Scientific  Officer  for  this  effort  and  his  excellent  work 
in  this  capacity  contributed  significantly  to  the  success  of  the 
program.  A  number  of  other  individuals  and  their 

organizations  played  an  important  part  in  the  success  of  this 
investigation,  and  it  is  important  to  the  authors  that  these 
people  be  recognized.  In  alphabetical  order,  they  are: 

Dr.  Donald  C.  Daniel  -  AFATL 
Mr.  Ramond  A.  Deep  -  Army  MICOM 
Mr.  Vernon  0.  Hoehne  -  AFWAL 
Mr.  Dale  E.  Hutchins  -  NAVAIR 

Mr.  Charles  M.  Jackson,  Jr.  -  NASA/Langley  Research  Center 
Dr.  Lionel  Pasiuk  -  NAVSEA 

Mr.  Wallace  C.  Sawyer  -  NASA/Lanqlev  Research  Center 
Dr.  Leon  H.  Schindel  -  NSWC 
Mr.  Edward  Sears  -  AFATL 

Mr.  David  S.  Shaw  -  NASA/Langley  Research  Center 
Dr.  Donald  J.  Spring  -  Army  MICOM 
Mr.  William  C.  Volz  -  NAVAIR  retired 
Mr.  W.  D.  Washington  -  Army  MICOM 

Sincere  thanks  go  to  the  wind  tunnel  staff  at  both  the  NASA 
Ames  and  Langley  Research  Centers  for  their  cooperation  and  help 
during  the  test  program. 

A  special  acknowledgement  is  due  Dr.  Jack  N.  Nielsen,  Chief 
Scientist  at  the  NASA/Ames  Research  Center.  He  conceived  the 
idea  of  a  cooperative  program  involving  the  three  services  and 
NASA,  brought  together  the  key  individuals  necessary  to  make  it 
happen,  and  directed  the  work  in  its  early  stages  while  he  was 
President  of  NEAR,  Inc. 

It  is  also  noteworthy  that  Dr.  Michael  J.  Hemsch,  co-author, 
was  intimately  involved  with  the  succesful  conduct  of  the  work 
while  at  NEAR,  Inc.,  and,  even  after  moving  to  PRC  Kentron,  he 
played  an  essential  role  in  the  comoletion  of  the  effort. 

Michael  R.  Mendenhall 
Auaust  1,  1986 
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Under  Office  of  Naval  Research  Contract  No.  N00014-80-C- 
0700,  Nielsen  Engineering  &  Research,  Inc.,  (NEAR)  conducted  an 
investigation  to  obtain  an  extensive  missile-fin  data  base  for 
use  in  broadly  applicable  engineering  prediction  programs  for 
calculating  the  aerodynamic  characteristics  of  body-tail  and 
canard-body-tail  missiles.  The  data  base  was  incoroorated  into  a 
specific  computer  program  called  PROGRAM  MISSILE  3  which  is  valid 
for  angles  of  attack  up  to  45°,  arbitrary  roll  angles,  fin 
deflection  angles  between  -40°  and  40°,  Mach  numbers  between  0.6 
and  4.5,  and  fin  aspect  ratios  between  0.25  and  4.0. 
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The  investigation  reported  herein  is  a  tri-services  coopera¬ 
tive  effo.  t  sponsored  by  the  Air  Force,  Army,  and  Navy  with 
extensive  participation  by  NASA/Ames  Research  Center  and 
NASA/Langley  Research  Center.  The  first  year's  work  involved: 
( 1)  selection  of  the  test  model  design,  test  parameters  and  test¬ 
ing  sequence,  (2)  preliminary  investigation  of  the  optimum 
approach  for  data  handling,  that  is,  preparing  the  data  for  and 
incorporating  it  into  PROGRAM  MISSILE  3,  and  (3)  revising  the 
equivalent  angle-of-attack  formulation  (Refs.  1  and  2)  to  incor¬ 
porate  the  new  fin  deflection  data  base.  The  results  of  that 
work  are  described  in  the  first  year's  report  (Ref.  3).  The 
second  year's  work  consisted  of  (1)  support  of  the  ongoing  wind 
tunnel  tests,  (2)  preparation  for  processing  the  data  to  be 
incorporated  into  the  data  base,  (3)  continued  improvement  of  the 
methods  used  in  MISSILE  3,  and  (4)  continued  code  development. 
The  results  of  that  work  are  described  in  the  second  year's 
report  (Ref.  4).  The  third  year's  work  continued  the  activities 
of  the  second  year  and  are  contained  in  the  third  year  report 
(Ref.  5).  The  fourth  year's  work  consisted  of  organizing  the 
experimental  data  base,  implementing  a  revised  eauivalent  angle 
of  attack  formulation,  and  incorporating  both  of  these  into  the 
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computer  code  MISSILE  3.  This  report  summarizes  the  previous 
work  and  presents  the  new  code  MISSILE  3. 

An  engineering  data  base  method  for  performance  prediction 
and  preliminary  design  of-  missiles  with  cruciform  fin  sections 
has  been  developed.  The  method  uses  a  newly  available  systematic 
data  base  which  covers  a  Mach  number  range  from  0.6  to  4.5  and 
fin  aspect  ratios  from  0.25  to  4.0,  angles  of  attack  up  to  45°, 
and  arbitrary  roll  angles.  Et  employs  the  equivalent  angle  of 
attack  concept  to  include  the  effects  of  vorticity  and  geometric 
scaling.  The  eouivalent  ancle  of  attack  method  has  been  modified 
somewhat  from  that  used  in  the  previous  MISSILE  programs 
described  in  References  1  and  2. 

The  report  is  divided  into  two  major  parts:  Experimental 
Program  and  Technical  Approach,  and  a  User's  Manual.  The  Experi¬ 
mental  Program  and  Technical  Approach  sections  describe  the 
experimental  tests,  the  data  base,  the  equivalent  angle  of  attack 
methodology,  and  the  fin  and  body  vortex  models.  It  is  not  the 
purpose  of  this  report  to  present  and  document  the  entire  tri¬ 
service  data  base.  The  incorporation  of  the  data  base  into  the 
engineering  prediction  method  is  the  primary  objective.  The 
User's  Manual  section  provides  instruction  for  the  use  of  the 
FORTRAN  computer  code  MISLE3 .  This  includes  a  program  descrip¬ 
tion,  program  limitations,  error  messages,  input  preparation,  and 
sample  cases.  The  report  concludes  with  comparisons  to  indepen¬ 
dent  experimental  data,  conclusions  and  recommendations. 

OVERVIEW  OF  FXPERIMENTAL  PROGRAM 


The  overall  objective  of  the  test  program  was  to  obtain  a 
high  Quality  systematic  fin  force  and  moment  data  base  for  use  in 
hiqh  angle  of  attack  missile  aerodynamics  computer  programs  based 
on  engineering  methods.  The  data  obtained  fall  into  two  cate- 
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gories:  (1)  fin  loads  without  fin  deflection  for  all  fins  in  the 
data  base  (stability  data),  and  (2)  fin  loads  with  fin  deflection 
for  a  control  fin  set  (control  data).  The  model  design  and  test 
conditions  were  based  on  use  of  Langley  Remote  Control  Missile 
Roll  Rig  and  are  intended  to  reflect  as  completely  as  possible 
the  range  of  flight  conditions  and  fin  designs  employed  by  modern 
high-performance  missiles.  The  Langley  model  is  2.5  inches  in 
diameter  and  allows  independent,  remote  control  of  four  control 
surfaces  and  model  roll  angle.  By  interchanging  modular  compo¬ 
nents,  canard,  wing,  or  tail  control  may  be  used.  For  the  tests, 
a  0.25  inch  shell  was  added  to  the  body  to  accommodate  a  three- 
component  balance  and  associated  wiring  for  each  of  the  four 
fins.  Hence,  the  diameter  of  the  cylindrical  portion  of  the  body 
was  3.0  inches.  All  of  the  data  obtained  for  the  data  base  was 
for  a  body-tail  configuration. 

Multiple  tunnel  entries  were  made.  Intermediate  and  hiqh 
Mach  number  (2.5  -  4.5)  tests  were  conducted  in  the  NASA/Langley 
Unitary  Plan  Tunnel,  and  low  and  intermediate  Mach  number  (0.6  - 
2.0)  tests  were  conducted  in  the  6-  by  6-foot  Supersonic  Tunnel 
at  NASA/Ames.  These  tests  were  performed  between  1982  and  1984. 

GENERAL  DESCRIPTION  OF  MISSILE  3  CODE 


The  range  of  parameters  for  MISSILE  3  are  described  in  this 
section.  The  Mach'  number  range  of  the  data  base  is  from  0,6  to 
4.5.  Test  Mach  numbers  of  0.6,  0.8,  0.9,  and  1.2  provide 
detailed  loads  for  the  transonic  speed  regime.  This  greatly 
improves  the  predictions  in  the  transonic  regime  compared  to 
previous  MISSILE  programs  (Refs.  1  and  2)  which  contained  data 
at  Mro  =  0.8  and  1.2.  The  supersonic  Mach  numbers  in  the  data 
base  are:  1.5,  2.0,  2.5,  3.0,  3.5,  and  4.5. 
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The  body  angle  of  attack  range  is  0°  <  <  45°.  Test 
angles  of  attack  were:  a  =0/2,  5,  10,  ...,  40,  45°.  The  fin 
control  deflection  angles  varied  from  -40°  to  40°  in  ten-degree 
increments.  The  body  roll  angle  range  is  0°  to  90°.  The  data 
base  contains  loads  for  -90°  <  <j>^<  90°  in  ten-degree  increments, 
where  is  the  fin  orientation  angle. 

The  fin  aspect  ratio  range*  of  the  data  base  is  0.25  to  4; 
this  greatly  extends  the  range  allowed  by  previous  MISSILE 
codes.  The  test  fins  had  aspect  ratios  of  0.25,  0.5,  1,  2,  and 
4,  and  the  taper  ratio  range  is  from  0  to  1.  The  fins  with 
aspect  ratios  of  1,  2  and  4  were  deflectable.  Trailing  edge 
sweep  and/or  forward  swept  leading  edges  are  not  permitted.  The 
test  fins  had  body  radius  to  semispan  ratios  (a/sm)  of  0.5.  The 
method  and  the  data  base  consider  cruciform  sections  with  iden¬ 
tical  fins  only. 

Figure  1  shows  the  angle  of  attack,  aspect  ratio,  and  Mach 
number  range  of  the  systematic  data  base,  and  the  range  of 
parameters  allowed  by  MISSILE  3  are  summarized  below: 

0.6  <  M  <  4.5 

00 

0°  <  a  <  45° 
c 

0 0  <  4>  <  90° 

0.25  <  AR  <  4.0 

0.0  <  X  <  1.0 

4  fins  per  section 

all  fins  in  a  section  are  identical 

symmetrical  airfoil  sections  only 

no  fin  trailing  edge  sweep 

axisymmetric  bodies 

TABLE  1.  Range  of  Parameters 


*Fin  aspect  ratio  is  tRe  aspect  ratio  of  two  fins  joined  at  the 
rootchord. 


The  data  base  used  for  MISSILE  3  is  much  larger  than  that 
used  in  MISSILE  2A.  The  MISSILE  2A  data  base  contains  experi¬ 
mental  data  for  Mach  numbers  from  0.8  to  3.0  and  fin  aspect 

ratios  from  0.5  to  2.0  ( AR  to  3.5  for  Ha  <  1.2).  Descriptions 
and  development  of  the  MISSILE  series  analysis  methods  can  be 
found  in  References  1,  2,  6,  7,  and  8. 

The  eauivalent  angle-of-attack  concept  and  the  data  base 

used  in  MISSILE  3  predict  missile  performance  similar  to  MISSILE 
2A.  However,  the  range  of  applicability  is  much  larger,  so  per¬ 
formance  can  be  investigated  over  a  larger  operating  range. 

MISSILE  3  allows  multiple  Mach  number  cases  in  a  single  run,  a 
convenience  not  available  with  MISSILE  2A.  MISSILE  3  has  also 
proven  to  be  more  economical  than  MISSILE  2A,  even  though  a 

larger  data  base  is  included.  MISSILE  3  has  been  compared  to 
independent  experimental  data  for  several  geometries  and  oper¬ 
ating  conditions,  and  it  has  provided  satisfactory  results. 
Further  comparisons  with  experiment  are  needed  to  determine  the 
full  range  of  operation  of  the  method  and  to  locate  areas  which 
require  further  improvement. 

OVERVIEW  OF  THE  ENGINEERING  METHOD 

This  section  of  the  report  describes  the  data  base  and  meth¬ 
odology  used  by  MISSILE  3.  The  missile  configurations  allowed  by 
MISSILE  3  are  cruciform  canara-boay-ta il ,  canard-body,  and  body- 
tail  configurations.  The  missile  is  divided  into  four  sections 
for  modeling  purposes;  the  nose  and  forebody,  the  canard  section, 
the  afterbody,  and  the  tail  section  as  shown  in  Figure  2.  For  a 
body-tail  configuration,  everything  forward  of  the  fin  section  is 
the  forebody. 


The  coordinate  system  and  fin  numbering  system  used  by 
MISSILE  3  are  shown  in  Figure  3.  The  x-axis  is  along  the  body 
centerline  and  positive  aft,  and  the  z-axis  is  in  the  wind  plane 
and  normal  to  the  body  axis  (positive  up).  The  y-axis  is  perpen¬ 
dicular  to  the  plane  containing  the  body  axis  and  the  wind  vector 
and  is  positive  to  the  right  looking  forward.  The  fins  are  num¬ 
bered  clockwise,  viewed  from  the  rear,  with  fin  1  always  in  the 
first  quadrant;  this  is  opposite  to  previous  MISSILE  codes  which 
numbers  the  fins  counterclockwise.  The  normal  force  for  fins  1 
and  3  is  positive  up  and  to  the  left,  and  the  normal  force  for 
fins  2  and  4  is  positive  up  and  to  the  right.  (See  Figure 
3(c).)  Fin  deflections  are  positive  when  they  tend  to  increase 
the  fin  incidence,  and  hinge  moments  are  positive  when  they  tend 
to  rotate  the  fins  so  as  to  increase  the  fin  incidence.  Fin 
bending  moments  for  fins  1  and  2  are  positive  when  they  give 
positive  contributions  to  the  rolling  moment,  and  bendina  moments 
for  fins  3  and  4  are  positive  when  they  give  negative  contribu¬ 
tions  to  the  rollinq  moment. 

Data  Base  Description 

The  data  base  consists  of  fin  data  and  wing-alone  data.  The 
fin  data  contain  normal-force  coefficient,  axial  center-of- 
pressure,  and  spanwise  center-of-pressure  information  with  and 
without  fin  deflections  for  fins  with  aspect  ratios  of  0.25  to 
4.  The  fin  data  base  is  divided  into  stability  (no  control 
deflections)  and  control  data,  and  each  of  these  sections  are 
further  divided  into  normal-force  coefficient  and  center-of- 
pressure  locations.  The  wing-alone  data  contain  normal-force 
coefficient  information  as  a  function  of  angle  of  attack  and  Mach 
number.  The  fin  data  base  is  discussed  first,  followed  by  the 
wing-alone  data  base. 
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Fin  Stability  Data  Base 


The  data  for  the  fin  data  base  come  from  tests  performed  at 
NASA  Ames  and  NASA  Langley  Research  Centers  as  part  of  the  Tri- 
Services  effort  {Ref.  9).  The  forces  and  moments  for  the  fins 
depicted  in  Figure  4  were  obtained  by  force  balance  measure¬ 
ments.  All  test  fins  have  body  radius  to  fin  semi-span  ratios, 
a/sm,  of  0.5.  These  tests  provide  normal-force,  hinge-moment, 
and  bending-moment  coefficients  for  Mach  numbers  between  0.6  and 
4.5,  angles  of  attack  up  to  45°,  and  roll  angles  from  -90°  to 
90°.  The  hinge-moment  and  bending-moment  coefficients  were 
reduced  to  provide  the  fin  centers-of-pressures,  x/C^  and 
y/s.  Table  2  below  depicts  the  aspect  ratio  and  taper  ratio 
domain  of  the  stability  data  base. 


Aspect 

Ratio 


region  of 
interpolation 


TABLE  2.  Aspect  Ratio  and  Taper  Ratio  Range  of  the  Stability 
Data  Base 


The  numbers  in  Table  2  are  the  test  fin  designation  numbers  (Fiq. 
4(a)),  and  the  Roman  numerals  represent  regions  of  interpolation. 


The  test  Mach  numbers  are: 

Mm=  0.6,  O.fi,  0.9,  1.2,  1.5,  2.0,  2.5,  3.0,  3.5,  4.5. 
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The  test  angles  of  attack  are: 

ac  =  0,  2,  5/  10,  15,  20,  25,  30,  35,  40,  45°. 

The  fin  orientation  angles  are: 

$  ~  — 90,  -80 , . . . ,  — 10,  0,  10, •••,80,  9  0  0  • 

Fin  Normal-Force  Coefficient.-  The  fin  normal-force  coeffi¬ 
cient,  CNF,  is  a  function  of  aspect  ratio,  taper  ratio,  Mach 
number,  included  angle  of  attack,  and  angle  of  roll. 
Methods  were  developed  to  interpolate  in  these  parameters  to 
obtain  specific  values  of  CNB'(  AR ,  X,M  m,  a  ,  <{>) .  For  example, 

linear  interpolation  in  the  Mach  number  parameter  for  normal- 
force  coefficient  results  in 

r  i 

CN(MJ  =  CN  ( M  x |^CN(M2)  -  CN  ( M  j  )  I  (1) 

where  Mj  and  M2  are  the  Mach  numbers  in  the  data  base  which  en¬ 
close  the  desired  M  . 

00 

Interpolation  in  the  aspect  ratio  and  taper  ratio  directions 
is  considered  separately  in  the  six  distinct  regions  depicted  in 
Table  2,  and  the  general  formulas  for  interpolation  in  these 
regions  follow. 

In  Region  I,  the  value  of  CN  is  the  point  in  the  AR-X  domain 
which  lies  on  the  plane  formed  by  fins  12,  31,  and  32;  therefore, 

Reqion  I:  0.0  <  X  <  0.5;  0.25  <  AR  <  0.50 

CNj  =  CN32  -  2( 0. 5- X) ( CN32  -  CN31)  -  4( 0.5-AR) ( CN32-CN12)  (2) 
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In  Region  II,  the  value  of  CN  is  the  point  in  the  AR-X 
domain  which  lies  on  the  plane  formed  by  fins  12,  32,  and  33; 

there  fore , 

Region  II;  0.0  <  X  <  0.25;  0.25  <  AR  <  0.50 

CNit  =  CN32  +  2{ X-0.5) (CN33-CN32)  -  4 (  0 . 5- AR) ( CN32-CN12)  (3) 

In  Region  III,  the  value  of  CN  is  estimated  first  as  the 
point  in  the  AR-X  domain  which  lies  on  the  plane  formed  by  fins 
31,  32,  and  42,  and  secondly  as  the  point  in  the  AR-X  domain 
which  lies  on  the  plane  formed  by  fins  42,  52,  and  51.  Since  it 
is  unlikely  that  these  two  planes  intersect  at  the  desired  point, 
the  interpolated  value  of  CN  is  taken  as  an  AR-weighte>.  combina¬ 
tion  of  the  two  estimates.  This  method  was  chosen  ove.  a  three 
point  method,  as  described  above  for  Regions  I  and  II,  because  it 
uses  more  information  and  should  provide  a  better  result.  There¬ 
fore, 

Region  III;  0.0  <  X  <  0.5;  0.50  <  AR  <  2.0 

0_  AR 

CNin  =[  CN32  -  2(0. 5-X)  (CN32-CN31)  +  2(  AR-0. 5)  (  CN4  2-CN3  2)  )•—“ 

+  [  CN52  -  2(0.5  -X)  (CN52-CN51)  -  (  2-  AR)  (  CN52-CN4  2 ) 

1  •  O 


(4) 

In  Region  IV,  the  value  of  CN  is  estimated  first  as  the 
point  in  the  AR-X  domain  which  lies  on  the  plane  formed  by  fins 
32,  33,  and  42,  and  secondly  as  the  point  which  lies  on  the  plane 
formed  by  fins  42,  52,  and  53.  As  in  Region  III  the  interpolated 
value  of  CN  is  taken  as  an  AR-weighted  combination  of  the  two 
estimates.  Therefore, 


-9- 


■.wji ujcw.ti.ji-jii_-iJi.VKji  jk.vw_t  w.t  »  A  Kjt_.vv»->-tx.Tjv^  vt.Vjjy^A^i.v;VCL,v\rJV‘»jj'»<.LVl^>a/i<.vVl.>..>, 


Region  IV:  0.0  <  X  <  1.0;  0.50  <  AR  <  2.0 

9— AT? 

CNIV  =[ CN32  +  2  (  X-0.5)  (CN33-CN32)  +  2(  AR-0.5)  (CN42-CN32) 

+  f  CN52  +  2( X-0.5) (CN53-CN52)  -  (  2-AR)  ( CN52-CN42)  ]  — y--g^ 


(5) 

In  Region  V,  the  value  of  CN  is  the  point  in  the  AR-X  domain 
which  lies  on  the  plane  formed  by  fins  51,  52,  and  62.  Thus, 

Region  V:  0.0  <  X  <  0.5?  2.0  <  AR  <  4.0 

CNy  =  CN52  -  2(0. 5-X) (CN52-CN51)  +  0. 5 ( AR-2 ) ( CN62-CN52)  (6) 

In  Region  VI,  the  value  of  CN  is  the  point  in  the  AR-X 
domain  which  lies  on  the  plane  formed  by  fins  52,  53,  and  62. 

Thus, 

Region  VI:  0.0  <  X  <  1.0;  2.0  <  AR  <  4.0 

CNVI  =  CN52  +  2{  X-.5)  (CN53--CN52)  +  0. 5(  AP-2)  (  CN62-CN52)  (7) 

Figure  5  illustrates  representative  examples  of  the  normal- 
force  coefficient  versus  angle  of  attack  for  several  Mach  numbers 
and  roll  angles  for  fin  52.  Recause  of  the  overwhelming  size  of 
the  data  base,  only  limited  examples  are  presented  in  this 
report . 

Fin  Center-of-Pressure ♦-  To  predict  the  fin  hinge-  and 
bending-moment  coefficients  and  the  missile's  rolling-  and  pitch¬ 
ing-moment  coefficients,  it  is  necessary  to  determine  the  axial 
and  lateral  positions  of  the  fin  center-of-pressure.  When  body 
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vortex-induced  effects  on  the  loads  are  removed,  t'u2  fin  center- 
of-pressure  locations  depend  primarily  on  the  fin  normal-force 
coefficient  and  therefore  are  insensitive  to  roll  angle.  This  is 
discussed  below.  The  data  base  is  used  to  obtain  the  center-of- 
pressure  location  associated  with  nonvortex  loads.  The  effects 
of  vorticity  in  the  flow  field  are  added  with  the  equivalent 
angle  of  attack. 

The  test  results  provide  normal-force,  hinge-moment  and 
bending-moment  coefficients  for  each  test  fin  as  a  function  of 
M  ,  a  and  <)>.  These  data  are  reduced  to  obtain  x/C  and  y/s  in 

00  C  K 

the  following  manner: 


x(Mw,ac,«{>) 


[HL  _  CHM(H°°'  V*0 

cr"  cn(m-'V*) 


y(MM,c(c,<j>) 


C_„  (M  ,  a  ,  <j>) 
BM  “  c 

Cm(M  ,  a  ,  4) 
N  <*>'  c 


(8) 

(9) 


In  order  to  obtain  the  centers-of-pressures  as  a  function  of 
|cn|  without  vortex  effects,  fin  2  on  the  windward  side  of  the 
body  is  considered  because  the  vortex-induced  loads  are  quite 

small.  The  measured  center-of-pressure  for  this  case  is  almost 
entirely  free  of  vortex  effects.  A  piecewise  linear  least 
squares  fit  to  the  "x( M  „r  ac  t  <l> ) /Cp  and  y(M  „,  ac ,  ‘JO/s  versus  CN(Mro, 
data  is  used  to  obtain  a  single  curve  fit  of  the  fin  2  data 

for  <{>  =  0°  to  80°.  Figure  6  depicts  typical  fits  of  the  x/C  and 

y/s  data  points  which  are  composed  of  roll  angles  from  0°  to 

80°.  As  with  the  normal-force  coefficient,  only  limited  repre¬ 
sentative  examples  from  the  data  base  are  shown.  Generally  the 
center-of-pressure  correlated  very  well  with  (c^j  for  all  fins 
and  Mach  numbers. 
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With  the  dependence  on  ac  and  4>  reduced  to  a  dependence  on 
|cN|  ,  x/CR  and  y/s  depend  only  on  AR,  X,  Mm  and  |  | .  The  in¬ 
terpolation  procedure  in  the  directions  is  a  bicubic 

spline,  and  the  interpolation  procedure  in  the  aspect  ratio  and 
taper  ratio  directions  is  the  same  as  that  used  for  the  fin 
normal-force  coefficient  presented  in  the  previous  section. 

Fin  Control  Data  Base 

The  control  data  base  is  made  up  of  two  parts:  transonic 

norrnal-force  coefficient  data  and  the  effect  of  control  fin 

deflection  on  the  x  and  y  centers-of-pressures.  The  experimental 
tests  provided  normal-force,  hinge-moment,  and  bending-moment 
coefficients  for  Mach  numbers  between  0.6  and  4.5,  angles  of 

attack  up  to  45°,  roll  angles  from  -90°  to  90°,  and  deflection 

angles  of  -40°  to  40°.  The  hinge-moment  and  bending-moment  coef¬ 
ficients  provide  the  fin  centers-of-pressures,  x"/C  and  y/s. 

R 

Deflected  fin  data  was  only  obtained  for  fins  with  aspect  ratios 
greater  than  or  equal  to  1.0.  The  following  table  depicts 
the  aspect  ratio  and  taper  ratio  domain  of  the  control  data  base: 


Aspect 
Ratio , 


C  D  -  Region  of 
interpolation 


TABLE  3.  Aspect  Ratio  and  Taper  Ratio  Range  of  the  Control 
Data  Base 

The  numbers  in  Table  2  are  the  test  fin  designation  numbers. 
The  missing  fins  denoted  by  Roman  numerals  are  considered  with 
interpolation  schemes  to  be  described. 
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The  test  Mach  numbers  are: 


M  a  =  0.8/  1.2,  2.0,  3.0,  4.5. 


The  test  angles  of  attack  are: 
o  =  0,  2,  5,  10,  15,  20,  25,  30,  35,  40,  45°. 

The  fin  orientation  angles  are: 

=  -90,  "80,...,  "10,  0,  10,. ..,80,  90°. 

The  test  deflection  angles  are: 

6  =  -40,  -30,  ...,  30,  40°. 

Fin  Control  Transonic  Normal-Force  Coefficient.-  The  fin 
normal-force  coefficient,  CNF,  is  known  as  a  function  of  aspect 
ratio,  taper  ratio,  Mach  number,  included  angle  of  attack,  angle 
of  roll,  and  deflection  angle.  Methods  are  developed  to 
interpolate  in  these  directions  to  obtain  CNF(  AR,  X,M  w,  a  ,<{>,<S). 
Linear  interpolation  is  used  in  the  Mach  number  and  in  the 
deflection  angle  parameters.  Interpolation  in  the  aspect  ratio 
and  taper  ratio  directions  is  considered  in  four  regions  for  the 
control  data  as  depicted  in  Table  3,  and  the  general  formulas  for 
interpolation  in  these  regions  follow. 

In  Region  III,  the  value  of  CN  is  the  point  in  the  AR-  ?. 
domain  which  lies  on  the  plane  formed  by  fins  42,  51,  and  52; 
thus, 

Region  III:  0.0  <  X  <  0.5;  1.0  <  AR  <  2.0 

CNm  =  CN52  -  2(0.5— A)  (CN52-CN51)  -  (  2-AR)  ( CN52-CN42)  (10) 
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In  Region  IV/  the  value  of  CN  is  the  point  in  the  AR-X 

domain  which  lies  on  the  plane  formed  by  fins  42,  52,  and  53; 

therefore , 

Region  IV:  0.0  <  X  <  1.0;  1.0  <  AR  <2.0 

CNIV  =  CN52  +  2( X-0.5) (CN53-CN52)  -  ( 2-AR) ( CN52-CN42)  0.1) 

In  Region  V,  the  value  of  CN  is  the  point  in  the  AR-X  domain 
which  lies  on  the  plane  formed  by  fins  51,  52,  and  62;  therefore, 

Region  V:  0.0  <  X  <  0.5;  2.0  <  AR  <  4.0 

CNV  =  CN52  -  2(0. 5-X) (CN52-CN51)  +  0 . 5 ( AP-2 ) ( CN62-CN52)  (12) 

In  Region  VI,  the  value  of  CN  is  the  point  in  the  AR-X 

domain  which  lies  on  the  plane  formed  by  fins  52,  53,  and  62; 

thus, 

Region  VI:  0.0  <  X  <  1.0;  2.0  <  AR  <  4.0 

CNVI  =  CN52  +  2( X-0.5) (CN53-CN52)  +  0 . 5 ( AR-2 ) ( CN62-CN52)  (13) 

Figure  7  illustrates  representative  examples  of  the  fin 
normal-force  coefficient  versus  angle  of  attack  for  several  Mach 
numbers  and  deflection  angles. 

Fin  Control  Cen ter-o f-Pressure ♦-  Prediction  of  the  fin 
hinge-  and  bending-moment  coefficients  and  the  missile's  rolling- 
and  pitching-moment  coefficients  reguire  that  the  effect  of 
control  deflections  on  the  axial  and  lateral  position  of  the  fin 
center-o  f-pressure  be  determined.  The  data  base  is  used  to 
obtain  the  center-of-pressure  location  associated  with  nonvortex 
loads.  Ry  examining  fin  2  at  <j>  =  0,  the  effect  of  control 
deflection  on  the  center-of-pressure  locations  is  obtained. 
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The  test  results  provide  normal-force/  hinge-moment/  and 
bending-moment  coefficients  for  each  test  fin  as  a  function  of 
M„f  ac,  <j>,  and  6.  The  «J>  =  0  data  are  reduced  to  obtain  x*/CR 
and  y/s  in  the  following  manner: 


x(M 


vs> 


•R 


'HL 

C„ 


<WM-' 


VS) 


cntf 


(14) 


cbh(m-'“o'S) 

- • -  = 

The  cenfcet’-of-pressure  data  are  plotted  versus  |cn|  for  each 
fin.  The  -40°  <  6  <  40°  data  foe  each  fin  is  plotted  on  the  same 
graph  for  each  Mach  number,  and  these  data  are  fit  with  a  piece- 
wise  linear  curve.  Representative  examples  of  these  fits  are 
shown  in  Figure  8.  The  dependence  on  ac  and  6  is  reduced  to  a 

dependence  on  |cn|,  and  x/CR  and  y/s  depend  only  on  AR,  X, 
Mroand  j | .  The  interpolation  procedure  in  the  AR-X  parameters 
is  the  same  as  that  used  for  the  fin  transonic  normal-force 
coefficient  described  in  the  previous  section. 


Wing-Alone  Data  Base 

A  wing-alone  data  base,  reauired  by  the  eauivalent  angle  of 
attack  method,  is  composed  of  normal-force  coefficient  data  for  a 
range  of  angles  of  attack  and  Mach  number.  The  aspect  ratio  and 
taper  ratio  range  for  this  data  base  are  the  same  as  those  for 
the  fin  stability  data  base  described  above.  The  wing-alone 
curves  were  generated  from  various  sources  (Refs.  10  and  11)  and 
from  manipulation  of  the  tri-service  control  data  base  for  posi¬ 
tive  fin  deflection.  Representative  examples  of  the  wing-alone 
results  are  shown  in  Figure  9. 
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Nose  Vortex  Shedding 


The  forebody  vortex  model  used  in  the  previous  MISSILE  pre¬ 
diction  methods  is  described  in  Reference  2  and  is  used  by 
MISSILE  3.  It  is  a  semi-empirical  model  (Ref.  12)  which  is  based 

on  measurements  on  various  bodies  for  a  <  20°  and  M  <2.  Sev- 

c  00 

eral  investigators  have  demonstrated  that  the  body  vortices  for 
such  flow  are  symmetrical,  and  this  finding  is  incorporated  into 
the  low  angle  of  attack  model.  For  angles  of  attack  greater  than 
20°,  the  data  on  body  vortices  are  sparse;  consequently,  a 
method  based  on  a  multi-vortex  tracking  code  (Ref.  13)  was  used 
to  calculate  a  nose  vortex  numerical  data  base  for  bodies  at  high 
angles  of  attack.  These  results  are  included  in  MISSILE  2A  and 
they  are  currently  used  in  MISSILE  3.  Reference  1  describes  the 
development  of  the  nose  vortex  data  base. 

Afterbody  Vortex  Shedding  and  Tracking 

Following  the  general  vortex  cloud  approach  used  in  Refer¬ 
ences  13,  14,  and  15,  the  three-dimensional  steady  flow  over  the 
afterbody  is  assumed  to  be  equivalent  to  a  two-dimensional, 
unsteady,  separated  flow.  The  two-dimensional  solution  is 
carried  out  in  the  crossflow  plane  where  the  flow  about  the  body 
in  the  presence  of  discrete  vortices  is  obtained.  At  each  time 
steo,  corresponding  to  an  interval  of  the  length  on  the  body,  a 
new  vortex  pair  is  shed  into  the  flow  field  from  separation 
points.  The  discrete  vortices  forming  the  wake  are  allowed  to 
move  under  the  influences  of  the  freestream  flow,  the  body,  and 
the  other  vortices  (including  fin  vortices)  in  the  field. 

The  calculation  procedure  is  carried  out  in  the  following 
manner  along  the  afterbody  following  the  canard  section.  Start¬ 
ing  at  the  crossflow  plane  located  at  the  canard  trailing  edge, 
the  boundary  layer  in  the  crossflow  plane  is  forced  to  separate 


at  points  on  the  body  determined  from  the  crossflow  Mach 
number.  Figure  10  shows  the  variation  of  the  separation  angle 
with  crossflow  Mach  number  as  obtained  from  References  16,  17, 

and  18.  The  equations  for  the  separation  line  used  in  MISSILE  3 
are : 

if  M  sina  <  0.8  =  90°  (16) 

®  c  s 

if  .8  <  M  sina  <  1.2  8  =  20°  (2.5  H  sina  -2)  +  90°  (17) 

«o  c  S  to  c 

if  M  sina  >  1,2  9  =  110°  (18) 

"  c  s 

The  trailing  vortices  from  the  canard  fins  are  placed  in  the 
flow  field,  and  vortices  are  shed  from  the  afterbody  at  the  pre¬ 
dicted  separation  points.  The  strength  of  the  body  separation 
vortices  is  determined  from  the  vorticity  transport  in  the  bound¬ 
ary  layer.  The  paths  taken  by  these  free  vortices  are  calculated 
by  integration  of  the  eauations  of  motion  of  each  vortex  in  a 
stepwise  fashion  using  a  variable  step  size  differential  equation 
solver.  A  detailed  development  of  the  multiple  vortex  tracking 
procedure  can  be  found  in  References  13  and  14.  At  following 
downstream  crossflow  planes,  the  canard  vortices  and  the  free 
vortices  being  shed  from  the  body  are  allowed  to  influence  the 
body  pressure  distribution  and  the  motion  of  other  vortices  in 
the  field.  This  procedure  is  carried  out  over  the  entire  length 
of  the  afterbody;  to  the  tail  leading  edge  for  a  canard-body-tail 
configuration  and  to  the  base  for  a  canard-body  configuration. 

Since  this  method  is  only  used  on  the  afterbody  where  the 
body  is  required  to  be  a  constant  radius  cylinder,  the  effects  of 
changing  body  radius  or  noncyl indr ical  bodies  are  not  consid¬ 
ered.  These  effects  could  be  included  in  MISSILE  3  at  a  later 
date  if  it  is  desirable. 
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Equivalent  Angle  of  Attack 

Data  bases,  such  as  the  one  used  in  MISSILE  3,  provide  a 
foundation  for  preliminary  design  and  prediction  methods,  but 
they  do  not  encompass  all  geometries  and  flow  conditions.  For 
example,  all  the  fins  in  the  new  data  base  have  radius  to  semi¬ 
span  ratios,  a/sm,  of  0.5  and  they  are  located  at  various  posi¬ 
tions  on  the  test  body.  In  order  to  investigate  missiles  with 
different  radius  to  span  ratios  and  fin  locations  different  from 
those  in  the  data  base,  it  is  necessary  to  have  a  method  which 
geometrically  scales  the  data  base  and  takes  into  account  dif¬ 
ferent  vortical  flow  fields.  The  equivalent  angle  of  attack 
concept  is  used  in  MISSILE  3  to  accomplish  these  tasks,  and  it  is 
discussed  briefly  in  this  section.  The  equivalent  angle  of 
attack  concept  is  presented  in  References  19  and  20;  however  the 
current  method  is  slightly  different  from  the  previous  methods. 

Scaling  Effects 

The  first  step  in  the  equivalent  angle  of  attack  approach  is 
to  determine  the  fin  normal-force  coefficient  from  the  data  base 
for  each  AP,  X,  M  ,  a  ,  and  4>.  combination  of  interest.  This  is 

00  C  X 

done  for  each  of  the  four  fins  and  is  denoted  as  CNF;  _.  The 
subscript  "i"  denotes  fin  number,  i=l,2,3,4,  and  the  subscript 
"o"  indicates  the  value  corresponds  to  a/sm  =  0.5.  If  the  cross- 
flow  Mach  number,  M^sina  ,  exceeds  unity,  it  is  necessary  to 
adjust  the  reference  of  CNFj  to  t^ie  averaae  local  dynamic 
pressure,  o^  .  This  correction  is  necessary  because  of  the 
presence  of  1&°  body  bow  shock  which  significantly  changes  the 
dynamic  pressure  near  the  fins.  Details  are  presented  in  the 
Appendix.  The  correction  is  based  on  predicted  flow  field 
results  from  the  Fuler  code  SWINT  {Ref.  21). 


The  next  step  is  to  determine  the  equivalent  angle  of  attack 


for  the  fin  in  the  data  base,  a 


If  the  crossflow  Mach 


number  exceeds  unity,  it  is  necessary  to  use  the  fin  average 

local  Mach  number,  M^  ,  because  of  the  presence  of  the 

i,o 

bow  shock.  As  above,  the  correction  is  based  on  predicted  flow 
field  results  from  program  SWINT.  For  the  fin  local  Mach 
number  and  the  normal-force  coefficient  of  interest  given  by 

CNF-  »(q  /q  ) .  ,  the  equivalent  anale  of  attack,  a  ,  is 

i/°  <*  x.  1,0  ■  eq. 

1,0 

calculated  from  the  wing-alone  tables  as  indicated  in  the 
following  sketch. 


CNF.  •  — 

i,o  - 


*£  . 

1,0 


Once  a  is  determined,  the  next  step  is  to  remove  the 

eq . 

1,0 

vortex-induced  effects  for  the  test  conditions  of  the  data  base; 
that  is,  the  nose  and  forebody  flow  field  and  the  body  radius  to 


semi-span  ratio,  a/sm  •  Tnese  adjustments  are  necessary  so  that 
the  vortex  field  for  the  actual  geometry  and  flow  conditions  can 
be  included  later.  The  fins  in  the  data  base  are  located  on  a 
cylindrical  body  at  various  axial  stations  aft  of  a  3-caliber 
ogive  nose.  A  semi-empirical  method,  based  on  3-caliber  ogive 
nose  and  cylinder,  is  used  to  obtain  the  body  vortex  strength  and 
location  for  the  test  fin  locations.  This  method  is  the  same  as 
that  used  for  the  nose  vortex  model  described  earlier.  With  the 
strength  and  location  of  the  forebody  vortex  known,  the  effect  of 
the  vortex  field  on  the  equivalent  angle  of  attack  is  obtained  in 
the  same  manner  as  the  previous  MISSILE  programs.  The  increment 
in  the  equivalent  angle  of  attack  due  to  the  body  vortex  system 
is  obtained  from  reverse  flow  theorems  as  described  in  Reference 
1  (Appendix  B).  The  following  equation  is  solved  for  the  vortex 

A 

free  equivalent  angle  of  attack,  a  ,  for  each  fin  <}>.  . 

eq.  '  Yi 

I/O 

A 

tan  cc  =  tana  -  tan(Aa  )  (19) 

qi ,o  qi ,o  eq  v, i ,o 

Once  the  vortex-free  equivalent  angle  of  attack  is  found, 
the  results  are  scaled  to  the  a/sm  of  the  actual  fin.  If  it  is 
assumed  that  Kv/,  the  Beskin  upwash  factor,  is  given  by 


A 

tan  a _ 

eq . 

K  =  r - (  20  ) 

w  tan  a  cos  4> 
c 

and  that  Kw  is  linear  with  respect  to  a/sm  as  predicted  by 
slender  body  theory.  Then  the  scaling  with  a/sm  is  given  by 

'  a  a 

tana  =  2  —  tana  -  (2  —  -  ])tana  cos<|>.  (21) 

eq .  ,  s  ea .  s  c  l 

i,l  m  i,o  m 

The  subscript  1  denotes  the  a/sm  of  interest,  and  the  ~  indicates 

the  absence  of  vortices.  With  a  known,  the  actual  vortical 

ea .  . 


flow  field  may  be  included.  If  the  first  fin  section  is  under 
consideration,  the  vortex  field  due  to  the  nose  and  forefcody  must 
be  added.  If  the  second  fin  set  is  being  considered,  the  vortex 
field  from  both  the  first  fin  set  and  the  afterbody  must  be 
included.  The  equivalent  angle  of  attack  for  the  fin  geometry 
and  flow  conditions  of  interest  for  no  fin  deflections  is  given 
by 

tana  =  tana  +  tan(Aa  )  (22) 

eqi,l  eqi,l  ea  v,i,l 

The  next  step  is  to  determine  the  fin  loads  without  deflec¬ 
tions.  To  calculate  the  fin  normal  force,  the  wing-alone  normal- 

force  coefficient  is  found  corresponding  to  a  and  M 

®qi,l  i,l 

Quantity  M  is  the  average  local  Mach  number  for  the  a/sm 
i  / 1 

of  interest.  With  CNWifi  and  the  average  local  dynamic 
pressure,  q  ^  ,  the  fin  normal-force  coefficient  is  calculated 


CNFi,l  ■  CNWi,l 


(23) 


The  bending-moment  coefficient  acting  on  the  fins  with  no 
control  effects  is  calculated  from  the  normal-force  coefficient 
and  the  y/s  center  of  pressure  in  the  stability  data  base.  The 

A 

vortex-free  fin  normal-force  coefficient,  CNF.  ,,  is  calculated 

a.  1/1 

is  used.  From  the 


in 


cne  same 


ma 


nner  as  CNF*  except  a 
i  *  i  ( 


ea 


i  1 1 


stability  spanwise  center-of-pressure  data  base,  (y/s).  ,  is 

i,i 


found  corresponding  to  CNW^  ^  and 


i  r  1 


for  the  fin  geometry  of 


interest.  The  bending-moment  coefficient  is  given  by 


CBMF . 


i,l  =  CNFi,l(^s)i,l  +  (CNFi,l  -  CNFifl)(y/s) 


v,i 


(24) 


where  (y/s) 

1  v ,  1 

the  vorticity 
theorems  (Ref. 


is  the  point  of  action  of  the  normal  force  due  to 
in  the  flow  field  as  determined  by  reverse  flow 
1). 
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The  hinge-moment  coefficient  acting  on  the  fins  with  no 
control  effects  is  calculated  using  the  normal-force  coefficient 
and  the  )c/Cn  from  the  stability  data  base.  (jT/C  ).  is  found 

K  a  K  i. 

from  the  data  base  corresponding  to  CNW.  .  and  M  .  The 

1,1  i,  1 

effect  of  vortices  on  the  axial  center-of-pressure  is  found  by 
the  method  of  Nielsen  and  Goodwin  {Ref.  22).  This  method  assumes 
the  effects  of  vortices  on  the  center-of-pressure  cause  it  to 
move  along  a  lifting  line  (constant  percent  chord). 
Given  (y/s)  .  from  reverse  flow  theorems,  (x/C  )  .  is  given  by 

V  f  1  K  V  f  1 


where 


(x/c_)  =  T  +  ( 1-  t)  f  ]-  X)  (y/s)  . 

K  V  f  X  v  /  1 


(X/Cp)ifl  _  (l-X)(y/s)lfl 

)  -  ( 1- X) (y/s) .  . 


(25) 


The  hinae-moment  coefficients  for  the  fins  are  calculated  as 


CHMF . 


J1L  _  (2 L_)  CNF 

CR  ‘Cp'i,!  CNFi,l 


+  —  ~  (£-)„  ,  (CNF.  .  -  CNF.  )  (27) 

CR  lr  Vfl  1,1  1,1 


This  completes  the  equivalent  angle  of  attack  formulation  and  the 
corresponding  fin  loads  for  no  control  deflections. 


Fin  Deflection  Effects 

The  effect  of  fin  deflection  on  the  equivalent  angle  of 
attack  is  presented  in  this  section.  For  Mach  numbers  greater 
than  1.5,  interference  effects  and  control  effectiveness  are 
handled  by  slender  body  theory  and  proper  use  of  o and  M  For 
Mach  numbers  less  than  1.2,  control  effectiveness  is  obtained 
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from  the  transonic  normal-force  coefficient  control  data  base, 
and  interference  effects  between  fins  are  obtained  from  slender 
body  theory.  For  1.2  <Moo  <  1.5,  a  combination  of  the  two 
methods  is  used.  The  procedure  used  herein  is  justified  by  the 
correlations  presented  in  Reference  19. 

The  equivalent  angles  of  attack  with  fin  deflection  for  each 
fin  for  M m  >  1.5  are  given  by 


*eq  1 , 2  ~  °eql  1  +  AsELF<Sl  ~  AADJ(q  /q  2, 1  XADJ62 

+  AOPP(q/q”)3,lXOPP53  +  AADJ(q/q“)4, 1XADJ64 


aeq2,2  °eq2  1  ~  1  XAdJ61  +  ASELF62 


+  AADJ(q/q~)3,lXADJ63  +  AOPP(q/q»)4,lXOPP64 


(28) 


a  =  a 

eg.  _  eq. 


+  AOPP(a/q~)  1,  1  XOPP6l  +  AADJ(q/q»)2,lXADJ62 


+  ASELF63  “  AADJ|(a  /q  4, 1  AADJ64 


*eq^ 2  aeq4  1  +  AftDJ(q  ^/q“5  If  1  XADJ6]  +  AOPP(q<t/q«’)2,lXOPP62 

AADJ(q/a»)3,lXADJ63  +  ASF,LF54 


(  30) 


For  M  <1.2,  the  equivalent  angle  of  attack  due  to  control 
deflection  is  given  by  the  previous  expressions  except  the 
Ac„rp<5.  terms  are  replaced  by  Aa  terms  which  are  calculated 

oijIjC  1  60  •  « 

l ,  6 

from  the  transonic  normal-force  coefficient  control  data 


base.  To  calculate  Aa  ,  the  normal-force  coefficients  with 

60  •  jf 

if 5 

and  without  deflection,  CNF.  .  and  CNF.  x  are  found  from 

i/O,  i,o,  6^=0 f 

the  data  base  for  the  H  ,  a  ,  A_,  X,  and  6.  combination  of 

«  c  R  l 

interest.  These  values  are  corrected  to  the  average  local  test 

dynamic  pressure  as  previously  described.  The  corrected  values 

of  CNF  and  the  average  local  test  Mach  number  are  used  to 

calculate  the  equivalent  angles  of  attack,  a  and 

eqi,« 

a  .  The  change  in  the  equivalent  angle  of  attack  of  fin  i 

eqi,6=0 

due  to  its  deflection  is  scaled  with  a/sm  using  slender-body 


theory  and  is  given  by 


=  (a 


SELF, a/s 
) _ ■■  ■■!" . 

!i,6=0  ASELF,a/sm=  .5 


Given  the  equivalent  angle  of  attack,  the  loads  on  the  fins 

are  calculated.  The  normal-force  coefficient  is  calculated  in 

the  same  manner  as  the  nondeflected  normal-force  coefficient. 

The  wing-alone  normal-force  coefficient  is  found  corresponding 

to  a  and  the  average  fin  local  Mach  number,  M  .  With 

eqi,2  _  *i,l 

CNW  •  and  the  average  local  dynamic  pressure,  q  ,  the  fin 

*•  t  X. .  , 

normal-force  coefficient  is  ' 


\,2  -  ^i^V^i,] 


The  effect  of  control  deflection  on  the  bending-  and  hinge- 
moment  coefficients  is  calculated  from  the  control  y/s  and  x/CR 
data  bases.  The  wing-alone  vortex-free  normal-force  coefficient 
for  the  deflected  fin  is  approximated  by 


CNW .  .  =  CNW .  .  -  ( CNW .  .  -  CNW .  . ) 

1,2  1,2  1,1  1,1 
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From  the  control  spanwise  center-of-pressure  data  base, 

A 

(y/sk-  ?  _  is  found  corresponding  to  CNW.  0  and  M  for  the  fin 

*(  ‘(S'  If  ‘ 

of  interest.  ' 

<Y/s>.,2  =  t(y/.)if2fC  -  <y/*>1,lfCl  +  (?/■>!,!  <351 


where  (y/s).  .  is  the  control  spanwise  center  of  pressure  cor- 
1  ,  l ,  c 

A 

responding  to  CNW.  .  and  M  .  (y/s).  „  is  the  spanwise  fin 

1,1  x.  •  1  f  ^ 

1 f  1 

center  of  pressure  with  control  and  vortex  effects  included. 


The  bending-moment  coefficient  is 

CBMFi,2  -  <™Fi,2(7/s>i,2  -  +  CBMFi,l 


(36) 


where  from  Equation  (33) 


CNFi,2  = 


The  hinge-moment  coefficients  are  calculated  from  the  con- 

is  found  corresponding 

t  C 

,  is  found  corresponding 


trol  x/CD 

A  *> 

tables 

• 

<  X/CR> i 

to  CNW.  1  and 
i/ 1 

A 

*i.  / 

if  1 

and 

<x/CR) 

to  CNW.  0  and 

1  f  l 

M  » 

The 

axial 

l.  . 

if  1 

flection  is  given  by 


<*/Vi,2  ■  «^CR'i(2,c  -  +  (x/Vi,l  (37> 


The  hinge-moment  coefficient  is  then  given  by 


C™Fi,2  *  (XHL/CR  -  <X/CR>i,2!  CNFi,2 


(  38) 


-  lxHL/CR  -  ‘’‘'Vi,!1  ™Fi,l  +  CSMFi,l 


This  completes  the  equivalent  angle  of  attack  method  imple¬ 
mented  within  MISSILE  3. 

Body  Loads 

The  body  force  calculation  used  in  MISSILE  3  is  the  same  as 
that  used  in  MISSILE  2A.  These  calculations  provide  a  very  rough 
estimate  to  the  forces  and  moments,  and  it  is  recommended  the 
user  input  body  forces  if  better  estimates  are  available. 

The  nose  forces  on  the  body  are  found  from  a  combination  of 
potential  and  crossflow  drag  theories.  Potential  forces  are 
included  by  specifying  the  normal-force  coefficient  slope, 
dC^/da.  The  crossflow  drag  term  is  included  to  approximate  vis¬ 
cous  and  vortex  loads.  The  formulation  for  these  loads  is  found 
in  Reference  4. 

The  afterbody  forces  are  calculated  by  a  combination  of 
crossflow  drag  theory  and  pressure  integrations.  The  crossflow 
drag  term  is  used  to  approximate  the  normal  force,  and  the  side 
force  is  calculated  by  integration  of  the  pressure  distribution 
on  the  afterbody.  These  formulations  are  presented  in  Reference 
2. 

MISLE3  DESCRIPTION  AND  USER'S  GUIDE 

MISLE3  is  an  aerodynamic  engineering  prediction  method  for 
preliminary  design  of  axisymmetric  missiles  with  one  or  two  cru¬ 
ciform  fin  sections.  The  foundations  of  the  method  are  an  exper¬ 
imental  data  base  and  the  equivalent  angle  of  attack  concept. 
The  equivalent  angle  of  attack  methodology  manipulates  the  data 
base  results  to  predict  fin  loadings  for  configurations  and  flow 
conditions  not  in  the  data  base.  The  range  of  parameters  per¬ 
mitted  in  MISLE3  are:  0°  <  a  <  45°,  0.6  <  M  <  4.5, 

c  00 

-40°  <  5  <  40°,  and  0°  <  <S>  <  90°. 
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The  following  sections  are  a  user’s  manual  for  the  program 
MISLE3.  Included  are  descriptions  of  the  program,  the  calcula¬ 
tion  procedure,  input  preparation,  sample  input  and  output  cases 
which  illustrate  various  program  features,  and  program  limita¬ 
tions. 

Program  Description 

This  section  provides  a  general  description  of  program 
MISLE3  and  its  various  subroutines.  The  code  is  written  in  modu¬ 
lar  form  so  that  as  improvements  or  additions  to  the  data  base  or 
methodology  become  available  they  are  easy  to  incorporate.  The 
flow  of  the  calculation  is  described  herein,  and  detailed  com¬ 
ments  are  provided  throughout  the  code  to  assist  the  user  in 

understanding  the  order  of  calculation. 

The  computer  code  consists  of  the  program  M  ISLE  3  and  60 
subprograms.  The  overall  flow  map  of  the  program  is  shown  in 
Figure  11  where  the  general  relationship  between  the  subroutines 
and  external  references  can  be  seen.  Communication  between  the 
program  modules  is  handled  by  argument  lists  and  named  common 

blocks.  A  cross  reference  table  for  the  calling  relationship 
between  the  program  subroutines  is  shown  in  Figure  12,  and  a 
similar  table  for  the  named  common  blocks  is  shown  in  Figure  13. 

The  program  is  written  in  standard  FORTRAN  77.  Execution 
for  a  single  case  on  a  VAX  11/750  can  vary  from  1.5  minutes  to 
mere  than  5  minutes  depending  on  the  geometry,  the  flow  condi¬ 
tions,  and  the  number  of  shed  vortices  on  the  afterbody.  MISLE3 
spends  a  significant  amount  of  time  setting  up  tables  in  which  to 
interpolate.  This  time  is  independent  of  the  number  of  cases 

beina  run;  therefore,  the  code  is  much  more  efficient  for  mul- 


Subroutine  Description 


This  section  briefly  describes  the  main  program  MISLE3  and 

its  various  subroutines. 

MISLE3  The  main  program  for  calculating  the  forces  and 

moments  on  missiles  with  cruciform  fin  sections. 

ALFEO  Interpolates  in  the  wing-alone , normal-force  coeffi¬ 

cient  table  to  determine  the  equivalent  angle  of 
attack  corresponding  to  CNF  and  the  Mach  number  of 
interest. 

BLOCK  DATA  Initializes  the  named  common  blocks  CNFD,  CNFS  and 
PREPRO. 

BODY 3  Sheds  and  tracks  vortices  in  the  midbody  section  at 

high  angles  of  attack  and  computes  loads  on  the 
midbody  section. 

BVTEX  Computes  symmetric  body  nose  vortex  strength  and 

position. 

CCL  Computes  the  slender-body  theory  span  loading  for 

fin-fin  interference. 

CHRT8  Calculates  B(dCL/dcc)  for  supersonic  Mach  numbers? 

Chart  8  of  NACA  Report  1307,  Reference  23. 

CH14Y6  Calculates  the  body  center  of  pressure  due  to  a 

wing  or  tail;  Charts  14,  15  and  16  of  NACA  Report 

1307,  Reference  23. 


m 


i  1 
$  u 


I 


i 


«r* 

& 


i 


EQ4 


E05 


E06 


E07 


DFEOKM 


FINLDS 
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FTRACK 


Calculates  kb(W)  +  ^  (sm  “ 
Equation  4  of  AIAA  Journal,  Vol. 
1982,  Page  856,  Reference  24. 

Calculates  kb(W)  +  ^  (sm  ” 
Equation  5  of  AIAA  Journal,  Vol. 
1982,  Page  856,  Reference  24. 


Calculates 


K 


*'B( w)  '“m 

Equation  6  of  AIAA  Journal,  Vol. 

1982,  Page  856,  Reference  24. 


8(1  +  X)  (s 


1}  (dCr/det)  from 
20,  No.  6,  June 


1)  (dCjVd a)  from 
20,  No.  6,  June 


1)  (dC^/d  a)  from 
20,  No.  6,  June 


Calculates 


K 


B(W) 


8(1  +  X)  (s  -  1)  (dCr/da)  from 
m  L 


Equation  7  of  AIAA  Journal,  Vol.  20,  No.  6,  June 
1982,  Page  856,  Reference  24. 

Integrates  a  system  of  N  first  order  differential 
eauations  from  X  to  X  +  DX  with  error  control  by 
using  a  Kutta-Merson  integration  scheme.  This 
routine  is  used  for  the  vortex  tracking  on  the 
a  f terbody. 

Computes  the  loads  on  a  fin  section. 

Computes  the  strengths  and  locations  of  the  trail¬ 
ing  vortices  from  the  first  fin  set. 

Calculates  the  derivatives  needed  for  the  ordinary 
differential  equation  solver  DFEQKM  for  tracking 
the  shed  vorticity. 


GVFL 


Computes  the  velocity  induced  by  a  set  of  NV  free 
vortices  at  an  arbitrary  field  point. 
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INFLU 


INTFAC 


INPT 


INTDG1 


Determines  the  region  of  influence  of  one  fin  of  a 

cruciform  set  on  another  for  M  >  1. 

00 

Calculates  the  fin-body  interference  factors,  KB^W) 
and  %(3)  and  the  body  center  of  pressure,  XB. 

Reads  in  and  prints  out  all  input  data  for  MISLE3. 

Computes  the  transformed  geometry  for  an  interdigi- 
tated  tail  fin  section.  Used  only  when  the  body 
bank  angle  exceeds  45°. 


INTDG2  Transforms  loads  for  interdigitated  tails  back  into 

the  original  coordinate  system.  Used  only  when  the 
body  bank  angle  exceeds  45°. 


LAMIJ 


I 


Computes  the  slender-body  theory  fin  deflection 
factors  given  the  slender-body  theory  span  loading 
for  fin-fin  interference. 


LNTRP  Performs  a  1-dimensional  linear  interpolation. 


LIN2D  Performs  a  2-dimensional  linear  interpolation. 


NOSE 

PROSS 

QLML 


Calculates  the  loads  on  the  body  nose. 

Sets  up  tables  of  wina-alone  normal-force  coeffi¬ 
cient  and  fin  centers  of  pressures  for  the  geometry 
of  interest. 

Calculates  (q_/q  )  and  Cm„/M  )  as  a  function  of 

Z  m  z  00 

M  sin(a  )  for  each  fin.  M  sin(a  )  >  1.  This  is 

00  C  00  C 

based  on  calculations  performed  using  the  Euler 
code  SWINT,  Reference  21. 
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REVFLO  Computes  the  weighted  average  of  the  vortex-induced 

upwash  on  each  fin  and  its  point  of  action. 

SETCNS  Interpolates  in  the  fin  stability  data  ba.>e  to 

produce  a  table  of  CNF(B)  for  no  fin  deflection  for 
the  geometries  of  interest  and  for  the  Mach  number 
of  interest.  Also  sets  up  a  table  of  transonic 
normal-force  coefficient  with  fin  deflection  if 
M  <  1.2. 

co 

SIMPl  Computes  an  integral  to  a  given  tolerance  using 

Simpson's  rule  and  a  given  function  for  the  inte¬ 
grand  . 

SIMSON  Computes  an  integral  using  Simpson's  rule  and  a 

given  array  of  ordinate  values. 

TRACK  Tracks  vortices  from  X  to  X  +  DX. 

VELBOD  Computes  the  velocity  induced  on  a  body  by  a  set  of 

NV  free  vortices. 

VORADD  Computes  the  initial  positions  and  strengths  of  the 

shed  body  vorticity. 

WNGCND  Interpolates  in  the  CNDxx  tables  (transonic  control 

normal-force  coefficient  data  base)  in  the  aspect 
ratio  and  taper  ratio  directions. 

Transonic  control  normal-force  coefficient  data  for 
FIN  42,  AR  =  1.0,  X  =  0.5,  as  a  function  of  M  , 


CND42 


.|*L  f *  (il.MTJ  AS* 


CND51  Transonic  control  normal-force  coefficient  data  for 

FIN  51/  AR  =  2.0/  X  =  0.0,  as  a  function  of  M  / 

00* 

ac/  $/  and  6. 

CND52  Transonic  control  normal-force  coefficient  data  for 

FIN  52/  AR  =  2.0/  X  =  0.5/  as  a  function  of  M  f 

CO 

o  ,  j>,  and  6. 

CND53  Transonic  control  normal-force  coefficient  data  for 

FIN  53/  AR  =  2.0/  X  =  1.0,  as  a  function  of  Mro, 
a  ,  4>,  and  <5. 

CND62  Transonic  control  normal-force  coefficient  data  for 

FIN  62/  AR  =  4.0/  X  =  0.5/  as  a  function  of  M  , 

00 

ac?  <}>/•  and  6. 

WNGCNT  Interpolates  in  the  CNTxx  tables  (stability  normal- 

force  coefficient  data  base)  in  the  aspect  ratio  j 

and  taper  ratio  directions.  1 

CNT12  Stability  normal-force  coefficient  data  for  FIN  12,  ! 

AR  =  0.25,  X  =  0.5,  as  a  function  of  M  ,  a  and  4>. 

OO  C 

CNT31  Stability  normal-force  coefficient  data  for  FIN  31, 

AP  =  0  c  5  -  X  =  n,0,  a  fnnct’On  M  n  6 

”  C  | 

CNT32  Stability  normal-force  coefficient  data  for  FIN  32, 

AR  =  0.5,  X  =  0.5,  as  a  function  of  M  ,  a  and  4> . 

00  C 

CNT33  Stability  normal-force  coefficient  data  for  FIN  33,  ! 

AR  =  0.5,  X  =  1.0,  as  a  function  of  M  ,  a  and  <J>. 

05  c  ( 

1 

CNT42  Stability  normal-force  coefficient  data  for  FIN  42, 

AR  =  1.0,  X  =  0.5,  as  a  function  of  M  ,  a  and  <j>.  i 

“  c  { 

j 

5 

f 

f 
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CNT51 


Stability  normal-force  coefficient  data  for  FIN  51, 


AR  =  2.0,  X  =  0.0,  as  a  function  of  M  ,  a  end  <j>. 

00  C 

CNT52  Stability  normal-force  coefficient  data  for  FIN  52, 

AR  =  2.0,  X  =  0.5,  as  a  function  of  M  ,  a  and  4>. 

CNT53  Stability  normal-force  coefficient  data  for  FIN  53, 

"AR  =  2.0,  X  =  1.0,  as  a  function  of  M  ,  o  and  $. 

»  c 

CNT62  Stability  normal-force  coefficient  data  for  FIN  62, 

AR  =  4.0,  X  =  0.5,  as  a  function  of  M  ,  o  and  <j>. 

WNGCNW  Constructs  a  table  of  wing-alone  normal-force  coef¬ 

ficients  as  a  function  of  M  and  a  . 

oo  c 

XBAR  Computes  the  vortex-free  axial  location  of  the  fin 

center  of  pressure  (stability  data  base). 

X8ARC  Computes  the  vortex-free  axial  location  of  the  fin 

center  of  pressure  for  fin  deflection. 

XCPWB  Calculates  the  body  longitudinal  center-of-pressure 

location  in  the  presence  of  wings  or  tails. 

BAP ( 1  +  X) ( 1  +  1/M  )  >  4  and  M  >  ]. 

GO  OO 

YPAR  Computes  the  vortex-free  spanwise  location  of  the 

fin  center  of  pressure  (stability  data  base). 

Computes  the  vortex-free  spanwise  location  of  the 
fin  center  of  pressure  for  fin  deflection. 
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Program  Restrictions  and  Limitations 

The  major  restrictions  in  MISLE3  were  discussed  previously, 
but  these  are  summarized  in  this  section  for  completeness.  Addi¬ 
tional  limitations  and  suggestions  based  on  the  authors*  expe¬ 
rience  with  the  code  are  also  included. 

The  scope  of  the  data  base  imposes  the  following  restric¬ 
tions  on  the  flow  conditions  and  the  geometries  that  can  be  ana¬ 
lyzed  : 

0.6  <  M  <  4.5 
00 

0°  <  a  <  45° 
c 

0°  <  4)  <  90° 

0.25  <  AR  <  4.0 

0  <  X  <  1 

4  fins  per  section 

All  fins  in  a  section  are  identical 
Symmetric  airfoils  only 
Mo  trailing  edge  sweep 
Axisymmetric  bodies 

Figure  1  illustrates  the  angle  of  attack,  aspect  ratio,  and  Mach 
number  ranee  of  the  systematic  data  base. 

while  the  data  base  in  MISLF3  is  limited  to  the  above  re¬ 
strictions,  the  program  will  allow  certain  of  these  restrictions 
to  be  violated,  but  it  will  print  out  a  warning  message.  These 
messages  are  discussed  in  the  following  section.  It  is  recom¬ 
mended  that  the  user  stay  within  the  aforementioned  limits  for 
general  usage  of  the  code,  and  any  excursions  outside  these 
limits  should  be  made  cautiously. 
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Error  Messages  and  Stops 

The  code  MISLE3  has  numerous  internal  error  messages  of 
varying  degrees  of  importance.  The  severest  error  message  is  a 
FATAL  which  results  in  a  numbered  program  STOP.  A  WARNING  mes¬ 
sage  is  printed  out  if  a  condition  is  violated  but  is  not  con¬ 
sidered  severe  enough  to  stop  execution  of  the  code  and  prevent 
the  remainder  of  the  cases  from  being  run.  The  error  messages 
and  their  location  in  the  code  are  described,  and  suggestions  to 
eliminate  the  problems  are  presented. 

STOP  10  This  is  a  FATAL  error  which  occurs  if  the  number  of 

input  Mach  numbers,  NMACH,  exceeds  5  or  the  number 
of  total  cases,  NMACH*NCOND,  exceeds  100.  One  of 
the  two  following  messages  will  be  printed: 

"FATAL  10  ****  NUMBER  OF  MACH  NUMBERS  EXCEEDS  5 

"FATAL  10  ****  TOTAL  NUMBER  OF  CASES  ( NMACH* NCOND) 
EXCEEDS  100  ****'' 

The  user  should  check  the  input  and  resubmit  the 
job.  Subroutine  INPT. 

STOP  20  This  is  a  FATAL  error  which  occurs  if  the  tip  chord 

of  a  fin  is  larger  than  the  root  chord.  One  of  the 
two  following  messages  will  be  printed  along  with 
the  quantities  in  question: 

"FATAL  20  ****  TIP  CHORD  EXCEEDS  ROOT  CHORD  FOR  FIN 
SET  1  ****" 

"FATAL  20  ****  TIP  CHORD  EXCEEDS  ROOT  CHORD  FOR  FIN 
SET  2  ****" 

Subroutine  INPT. 

STOP  30  This  is  a  FATAL  error  which  occurs  when  a  fin  is 

not  fully  on  the  body  or  when  the  rear  fin  section 


overlaps  the  front  fin  section.  One  of  the  four 
following  messages  is  printed  along  with  the  quan¬ 
tities  in  question: 

"FATAL  30  ****  LEADING  EDGE  OF  FIN  SET  1  LESS  THAN 
XZERO  ****" 

"FATAL  30  ****  TRAILING  EDGE  OF  FIN  SET  1  EXCEEDS 
XBASE" 

"FATAL  30  ****  LEADING  EDGE  OF  FIN  SET  2  LESS  THAN 
TRAILING  EDGE  OF  OF  FIN  SET  1  ****" 

"FATAL  30  ****  TRAILING  EDGE  OF  FIN  SET  2  EXCEEDS 
XBASE" 

Subroutine  INPT. 

STOP  40  This  FATAL  error  occurs  if  any  of  the  input  Mach 

numbers  are  less  than  0.6  or  greater  than  4.5.  One 
of  two  following  messages  is  printed: 

"FATAL  40  ****  MACH  NUMBER  LESS  THAN  .6  ****'' 

"FATAL  40  ****  MACH  NUMBER  EXCEEDS  4.5  ****" 

STOP  50  This  FATAL  error  occurs  if  any  of  the  input  fins 

have  aspect  ratios  less  than  0.25  or  exceeding 
4.0.  One  of  the  four  following  messages  is  printed 
along  with  the  quantities  in  question: 

"FATAL  50  ****  ASPECT  RATIO  FOR  FIN  SET  1  LESS  THAN 
.25  ****" 

"FATAL  50  ****  ASPECT  RATIO  FOR  FIN  SET  1  EXCEEDS 
4.0  ****•' 

"FATAL  50  ****  ASPECT  RATIO  FOR  FIN  SET  2  LESS  THAN 
.25  ****" 

"FATAL  50  ****  ASPECT  RATIO  FOR  FIN  SET  2  EXCEEDS 
4.0  *  *  *  * l! 

Program  MISLE3. 
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WARNING  60  This  WARNING  occurs  if  any  of  the  input  angles  of 
attack  exceed  45  degrees.  A  program  STOP  is  not 
used  here  since  the  user  may  be  running  a  series  of 
angles  of  attack.  The  program  will  skip  over  any 
conditions  for  which  this  warning  occurs.  The 
following  error  message  is  printed: 

"WARNING  60  ****  "ANGLE  OF  ATTACK  EXCEEDS  45 

DEGREES" 

Program  MISLE3 . 

WARNING  70  This  WARNING  occurs  if  any  of  the  input  angles  of 
attack  exceed  30  degrees  when  the  Mach  number  is 
less  than  0.8.  A  program  STOP  is  not  used  here 
since  the  user  may  be  running  a  series  of  angles  of 
attack  and  several  Mach  numbers.  The  program  will 
skip  over  conditions  for  which  this  warning 
occurs.  The  following  error  message  is  printed: 
"WARNING  40  ****  ANGLE  OF  ATTACK  EXCEEDS  30  DEGREES 
FOR  MACH  NUMBER  LESS  THAN  .8  ****".  The  points 

that  violate  the  conditions  are  printed. 

Program  MISLE3. 

WARNING  80  This  WARNING  occurs  if  the  aspect  ratio,  Mach  num¬ 
ber,  and  angle  of  attack  combination  are  outside 
the  data  base  values  illustrated  in  the  preceding 
sketch.  One  of  the  two  following  messages  is 

printed  along  with  the  parameters  of  interest: 
WARNING  80  ****  ASPECT  RATIO,  MACH  NUMBER,  ANGLE  OF 
ATTACK  COMBINATION  OUTSIDE  OF  DATA  BASE  FOR  FIN  SET 
1  *  *  *  *  « 

WARNING  80  ****  ASPECT  RATIO,  MACH  NUMBER,  ANGLE  OF 
ATTACK  COMBINATION  OUTSIDE  OF  DATA  BASE  FOR  FIN  SET 

2  *  *  *  *  it 

Program  MISLE3. 
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Input  Description 


This  section  describes  the  input  da'.e*  and  format  required  by 
program  MISLE3.  MISLE3  uses  list-directed  input  to  eliminate  the 
tediousness  of  input  preparation.  MISLF3  does  not  require  any 
specific  knowledge  of  the  dimensional  unit  of  measure,  but  the 
user  must  be  consistent  throughout  and  not  mix  units.  The 
remainder  of  this  section  describes  the  order  of  input  and  the 
input  variables. 

Item  1  is  a  single  card  containing  an  integer  specifying  the 
number  of  title  cards  to  be  input. 

Item  Variable  Description 

1  NCARD  Number  of  title  cards  used  to  identify 

the  run;  NCARD  >  1. 

Item  2  is  a  set  of  NCARD  title  cards  used  to  identify  the 
run.  These  cards  are  limited  to  an  80  column  field. 

Item  3  is  a  single  card  containing  the  program  control  vari¬ 
ables. 

Item  Variable  Description 

3  NFIN  Number  of  fin  sets;  NFIN  =  1  or  2 . 

NCOND  Number  of  attitude  conditions  (acr  ,  <$.. 

combinations)  for  which  calculations  are 
to  be  made;  NCOND  <  100.  Note  that 
NMACH*NCOND  <  100. 
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NMACH  Number  or  Mach  numbers  for  which  NCOND 

conditions  are  to  be  run;  NMACH  <  5.  Note 
that  NMACH*NCOND  <  100. 

NXAB  Number  of  axial  body  stations,  aft  of  the 

first  set  of  fins,  at  which  body  vortices 
are  to  be  shed.  Also,  the  number  of 
items  at  which  vortex  information  will  be 
output  if  0UTP=.T.  NXAB  <  40.  Note  that 
( 2*NXAB+4*NSHED)  <  100.  NXAB  should  be 
chosen  such  that  the  spacing  between 
stations  is  less  than  D/2. 

NSET  Number  of  the  fin  set  used  for  control. 

=  1,  First  fin  set  used  for  control 
=  2,  Second  fin  set  used  for  control 

NSHED  Number  of  trailing  vortices  shed  per  fin 

from  the  forward  (first)  set  of  fins; 

1  <  NSHED  <  10.  This  option  is  provided 
for  closely  coupled  canard  and  tail 
fins.  NSHED  =  1  is  sufficient  if  there 
is  a  long  afterbody  between  the  fin  sec¬ 
tions. 

LTAIL  Logical  variable  concering  interdigita- 

tion  of  the  rear  (second)  set  of  fins 
with  respect  to  the  first  set. 

=  F,  Second  fin  set  in  line  with  the 
first  set. 

=  T,  Second  fin  set  is  rolled  45°  clock¬ 
wise  (viewed  from  rear)  with  re¬ 
spect  to  the  first  set. 
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E5  Accuracy  criterion  for  vortex  tracking. 

Suggested  value,  E5  <  .001. 

BSHED  Logical  variable  concerning  shedding  of 

body  vortices  in  segment  after  the  first 
fin  set. 

=  F,  No  shedding  of  body  vortices. 

=  T,  Body  vortices  shed  at  the  NXAB 

stations  if  separation  conditions 
are  met. 

OUTP  Logical  variable  concerning  intermediate 

output.  If  OUTP  is  true  vortex  tracking 
information  for  diagnostic  purposes  is 
printed. 

=  F,  No  intermediate  output. 

=  T,  Intermediate  output  will  be 
printed . 

DEFLEC  Logical  variable  concerning  fin  deflec¬ 

tion  . 

=  F,  Fins  are  not  deflected. 

=  T ,  NSET  fin  set  has  deflected  fins. 

Item  4  is  a  single  card  containing  the  reference  information 
for  the  runs.  There  is  not  a  specific  unit  of  measure  used  in 
MISLE3,  but  the  user  must  be  consistent. 

Item  Variable  Description 

4  SROUT  Reference  area,  dimensional. 

LROUT  Reference  length,  dimensional. 

XMC  Moment  center  of  missile,  dimensional. 
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s 


s 

v> 


1 


(l! 


D 


XZERO 


XBASE 


Diameter  of  missile,  dimensional. 

Axial  location  of  the  nose  tip,  dimen¬ 
sional. 

Axial  location  of  missile  base,  dimen¬ 
sional. 


Item  5  is  a  single  card  containing  the  geometric  information 


for 

Si 

W 

the  first  set 

of  fins. 

Item 

or 

Variable 

Description 

h  5 

SP  AN(  1) 

Exposed  fin  semispan,  dimensional. 

1 

XLE( 1) 

Axial  location  of  the  leading  edge  of  the 

Lx« 

fin  root  chord,  dimensional. 

:  R 

HL(  1) 

Axial  distance  from  the  leading  edge  of 
the  fin  root  chord  to  the  fin  hinge  line. 

i 

dimensional . 

fl 

CT(  1) 

Fin  tip  chord,  dimensional. 

y 

CR(  1) 

Fin  root  chord,  dimensional. 

Item  6  is  a 

single  card  containing  the  geometric  information 

fe  for 

the  second  set  of  fins  and  is  included  only  if  NFIN  =  2. 

1 

Omit  item  6 

if  NFIN  £  2. 

ft  Item 

ft 

Variable 

Description 

^  6 

SPAN ( 2) 

Exposed  fin  semispan,  dimensional. 

g 
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XLE(2)  Axial  location  of  the  leading  edge  of  the 

fin  root  chord,  dimensional. 

HL(2)  Axial  distance  from  the  leading  edge  of 

the  fin  root  chord  to  the  fin  hinge  line, 
dimensional . 

Fin  tip  chord,  dimensional. 

Fin  root  chord,  dimensional. 

Item  7  is  a  set  of  NCOND  cards  specifying  the  attitude  con¬ 
ditions.  1  <  I  <  NCOND. 

Item  Variable  Description 

7  ALFAC(I)  Body  angle  of  attack;  angle  between  the 

wind  vector  and  the  body  axis  in  degrees. 
0°  <  a  <  45°. 

PHI ( I )  Bank  angle  in  degrees;  angle  between  the 

Z-axis  and  the  hinge  line  of  fin  1  of  the 
first  fin  set;  angle  measured  clockwise 
(looking  forward).  0°  <  <{>  <  90°. 

DELTA( 1,1)  Deflection  of  fin  1  in  degrees;  positive 

if  it  produces  a  counterclockwise  rolling 
moment  (looking  forward). 

-40°  <  6  <  40°. 

DF.LTA(2,I)  Deflection  of  fin  2  in  degrees;  positive 

if  it  produces  a  counterclockwise  rolling 
moment  (looking  forward). 


CT(  2) 
CR(  2) 


a&3M,sma^VKJ5eaa5m2a533cssssssm».\Kjci2^ 


DELTA ( 3/1)  Deflection  of  fin  3  in  degrees?  positive 
if  it  produces  a  clockwise  rolling  moment 
(looking  forward).  -40°  <  6  <  40°. 

DELTA(4,I)  Deflection  of  fin  4  in  degrees;  positive 
if  it  produces  a  clockwise  rolling  moment 
(looking  forward).  -40°  <  6  <  40°. 

Item  8  is  a  single  card  containing  the  Mach  numbers  for  the 
run;  1  <  I  <  NMACH. 


Item  Variable  Description 

8  FMACH(I)  Freestream  Mach  number;  1  <  I  <  NMACH 

0.6  <  M  <  4.5. 

00 

Item  9  is  a  single  card  containing  an  integer  flag,  I FLB , 
which  determines  whether  the  body  loads  are  input  or  calculated. 

Item  Variable  Description 

9  IFLB  =  0  Body  loads  are  input  in  item  10. 


=  1  Body  loads  are  calculated;  nose 

information  must  be  input  in  Items 
11  -  15. 


Item  10  is  a  set  of  NMACH*NCOND  cards  containing  the  input 
loads?  1  <  IJ  <  NMACH*NCOND,  one  set  of  loads  per  card.  These 
values  should  be  nondimensional ized  with  respect  to  the  input 
reference  area  and  length. 

Omit  item  10  if  IFLB  >  0. 
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Item 


Variable 


Description 


CZB(IJ) 


Z-direction  force  coefficient  of  the  body 
alone  {normal  force)  for  condition  IJf 
where  IJ  =  I+( J-l)*NCOND, 

1  <  I  <  NCOND,  1  <  J  <  NMACH.  The  inner 
loop  is  on  I  (attitude  conditions),  and 
the  outer  loop  is  on  J  (Mach  numbers). 


CMYB(IJ) 


Moment  coetficient  around  y-axis  of  the 
body  alone  (pitching  moment)  for  condi¬ 
tion  IJ. 


CXB(IJ) 


X-direction  force  coefficient  of  the  body 
alone  (axial  force)  for  condition  IJ. 


Omit  items  11-15  if  IFLB  <  0. 

Item  11  is  a  single  card  containing  an  integer  flag,  IDCN, 
which  determines  whether  the  potential  normal-force  coefficient 
slope  dCN/da  for  the  body  alone  is  input.  If  DCNDA  is  not  to  be 
input  in  Item  12,  it  is  set  to  its  default  value  of  2.0. 


Item 


Variable 


Description 


IDCN 


=  0  The  default  value  of  DCNDA  is  used, 
DCNDA  =  2.0. 


=  1  DCNDA  is  input  in  Item  12. 

Item  12  is  a  single  card  containing  the  potential  normal- 
force  coefficient  slope  for  the  body  alone,  DCNDA. 

Omit  Item  12  if  IDCN  <  0. 


Item 


Variable 


Description 


DCNDA 


Potential  normal-force  coefficient  slope, 
dCN/da,  in  radians. 


Item  13  is  a  single  card  containing  the  number  of  points 
specifying  the  nose  geometry,  NNOSE. 


Item 


Variable 


Description 


NNOSE 


Number  of  points  specifying  the  body 


nose. 


Items  14  and  15  are  two  sets  of  NNOSE  values  specifying  the 
body  nose;  XNOSE(I),  RNOSE(I),  1  <  I  <  NNOSE. 


Item 


Variable 


Description 


XNOSE(I) 


RNOSE( I) 


x-location  at  which  the  body  radius 
is  specified,  dimensional. 

1  <  I  <  NNOSE. 

Body  radius  at  XNOSE(I),  dimen¬ 
sional.  1  <  I  <  NNOSE. 


Sample  Input  Cases 


In  this  section,  two  sample  input  cases  are  described  to 
illustrate  the  various  proqram  options  available  in  MISLE3. 
Their  purpose  is  to  illustrate  the  available  futures  and  options 
and  to  provide  sample  inputs  which  will  help  users  prepare  input 
for  specific  cases.  An  effort  has  been  made  to  simplify  the 
input  requirement  of  MISSILE  3  compared  to  MISSILE  2 A.  Figure 
14  shows  the  aeometric  inputs  required  by  MISSILE  3. 
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Sample  Case  1,  Figure  15/  is  a  test  case  based  on  the  AIM9 
wind  tunnel  model  with  king  size  tail  fins  (Ref.  25).  This  case 
will  also  be  used  as  the  sample  output  case  in  the  next  sec¬ 
tion.  It  is  recommended  that  this  case  be  run  initially  to 
provide  a  check  with  the  results  presented  in  the  following  sec¬ 
tion. 

Input  Item  1/  NCARD,  indicates  two  title  cards  are  to  be 
input.  The  following  two  cards  are  the  reauired  title  cards  for 
the  run. 


Input  Item  3  contains  various  program  control  options.  NFIN 
=  2  indicates  that  the  missile  has  two  cruciform  fin  sections. 
NCOND  =  3  indicates  that  3  attitude  conditions  are  specified  in 
Item  7.  Results  for  these  three  conditions  are  calculated  for 
each  input  Mach  number.  NMACH  =  1  specifies  one  Mach  number  is 
input  in  Item  8.  Vortices  after  the  first  fin  set  are  shed  and 
tracked  at  25  stations/  NXAB,  along  the  afterbody  ahead  of  the 
second  fin  section.  NSET  =  ]  indicates  that  if  a  fin  set  is 
deflected  for  control  effects/  it  will  be  fin  set  1.  NSHED  =  1 
specifies  one  trailing  vortex  per  fin  is  shed  from  the  first  fin 
set.  The  second  fin  set  is  not  interdigitated  45  degrees  with 
respect  to  the  front  set,  LTAIL  =  F.  The  accuracy  criterion  for 
the  vortex  tracking,  E5,  is  0.001.  The  logical  variable  BSHED  is 
set  to  T  so  that  vortices  are  shed  from  the  afterbody  if  the 
separation  conditions  are  met.  OUTP  =  F  indicates  that  vortex 
tracking  output  along  the  afterbody  is  not  output.  The  final 
variable  in  Item  3,  DEFLEC  =  F,  indicates  there  are  no  deflected 
fins. 

Input  Item  4  contains  the  reference  information  for  the 
run.  These  items  are  dimensional,  and  for  this  case  the  unit  of 
measure  is  centimeters.  The  reference  area,  SROUT,  is  the  maxi¬ 
mum  crossectional  area  of  the  missile?  the  reference  length, 
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LROUT,  is  the  maximum  diameter;  and  the  moment  center,  XMC,  is 
46.04  centimeters  from  the  nose.  The  maximum  diameter  is  4.234 
centimeters.  The  axial  location  of  the  nose  tip  is  0.0,  and  the 
axial  location  of  the  missile  base  is  102.319  cm. 

Input  Item  5  contains  the  geometry  and  location  of  the  first 
fin  set.  The  exposed  semispan,  SPAN(l),  is  7.24  centimeters. 
The  root  chord  leading  edge,  XLE(l),  is  located  at  axial  station 
15.80,  and  HL(1)  indicates  the  hinge  line  is  6.35  centimeters  aft 
of  the  root  chord  leading  edge.  The  tip  and  root  chords,  CT( 1) 
and  CR( ] ) ,  are  3.338  and  11.123  centimeters,  respectively. 

Input  Item  6  is  the  same  as  item  five  except  the  geome¬ 
try  and  location  of  the  second  fin  set  is  input:  SPAN(2)  =  9.04, 
XLE( 2 )  =  80.729,  HL( 2)  =  10.795,  CT( 2)  --  12.55,  and  CR(2)  = 
12.59.  All  these  values  are  specified  in  centimeters.  The  hinge 
line  for  this  fin  set  is  arbitrary  because  it  will  not  be  de¬ 
flected  in  this  run. 

Input  Item  7  contains  the  NCOND  attitude  conditions  for 
which  results  are  to  be  obtained.  For  this  sample  case,  an  angle 
of  attack  sweep  is  demonstrated.  The  first  attitude  condition 
for  the  body  is  2.0°  angle  of  attack  and  45°  roll  angle.  There 
are  no  deflections  for  this  case,  but  values  must  be  input 
because  tno  code  uses  list-directed  input  and  expects  input 
values.  The  four  values  of  deflection  angles  are  input  0.0.  The 
second  c'.nd  third  attitude  conditions  are  identical  to  the  first 
except  the  angle  of  attack  is  changed  to  5.0°  and  10.0°,  respec¬ 
tively. 

Input  Item  8  contains  the  run  MACH  numbers.  For  this  case, 
a  single  Mach  nun.ber  is  set  to  2.50. 
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Input  Item  9,  IFLB,  is  a  single  integer  flag  indicating 
whether  the  body  forces  are  to  be  input  or  calculated  by 
MISLE3.  For  this  case  they  are  calculated,  IFLB  =  1.  Input  Item 
10  is  not  reauired  since  the  body  forces  are  calculated. 


Input  Item  11,  IDCN,  is  a  single  integer  flag  indicating 
whether  the  potential  normal-force  coefficient  slope,  dC^/da,  is 
input  or  defaulted  by  the  program.  IDCN  =  1  indicates  that  DCNDA 
is  to  be  input  in  Item  12.  Item  12,  DCNDA,  is  the  potential 
normal-force  coefficient  and  is  set  to  2.0. 

Items  13  through  15  contain  the  body  nose  geometry  required 
for  the  body  force  calculation.  NNOSE  =  31  in  Item  13  specifies 
31  points  are  used  to  define  the  body  nose.  The  31  axial  sta¬ 
tions,  XNOSE(I),  where  the  nose  radius  is  defined  are  input  in 

Item  14.  The  nose  radii,  RNOSE(I),  at  the  31  axial  station 
spejified  in  Item  14,  are  input  in  Item  15.  Item  15  completes 
the  input  for  Sample  Case  1. 

Sample  Case  2,  Figure  16,  is  a  test  case  based  on  the  Army 

Generalized  Missile  (Ref.  26)  for  multiple  Mach  numbers  and 

deflections.  Input  Item  ],  NCARD,  indicates  three  title  cards 
are  input.  The  following  three  cards  in  item  2  are  the  title 

cards  for  the  run. 


Input  Item  3  contains  various  program  control  options.  NFIN 
=  2  indicates  that  the  missile  has  two  cruciform  fin  sections. 
NCOND  =  4  specifies  four  attitude  conditions  are  input  in  item 
7.  Results  for  these  four  conditions  are  calculated  for  each 
input  Mach  number.  NMACH  =  2  indicates  two  Mach  numbers  are 
I  input  in  Item  8.  Body  vortices  after  the  first  fin  set  are  shed 

and  tracked  at  10  stations,  NXAB,  along  the  afterbody  ahead  of 
the  second  fin  section.  NSET  =  1  identifies  fin  set  1  as  control 
i  fins.  NSHED  =  1  indicates  one  trailing  vortex  per  fin  is  shed 


from  the  first  fin  set.  The  second  fin  set  is  not  interdigitated 
45  degrees  with  respect  to  the  front  set,  since  LTAIL  =  F.  The 
accuracy  criterion  for  the  vortex  tracking,  E5,  is  0.001.  The 
logical  variable  BSHED  is  true  so  that  vortices  are  shed  from  the 
afterbody  if  the  separation  conditions  are  met.  OUTP  =  T  indi¬ 
cates  the  vortex  tracking  input  along  the  afterbody  is  printed. 
The  final  variable  in  Item  3,  DEFLEC  =  T,  indicates  that  fins  of 
fin  set  1  are  deflected  for  control  effects. 

Input  Item  4  contains  the  reference  information  for  the 
run.  These  items  are  dimensional,  and  for  this  case  the  unit  of 
measure  is  inches.  The  reference  area,  SROUT,  is  the  maximum 
crossectional  area  of  the  missile,  the  reference  length,  LROUT, 
is  the  maximum  diameter,  and  the  moment  center,  XMC,  is  26.0 
inches  from  the  nose.  The  maximum  diameter  is  D  =  5.0  inches. 
The  axial  location  of  the  nose  tip  is  0.0,  and  the  axial  location 
of  the  missile  base  is  52.00  inches. 

Input  Item  5  contains  the  geometry  and  location  of  the  first 
fin  set.  The  exposed  semispan,  SPAN(l),  is  3.75  inches.  The 
root  chord  leading  edge,  XLE(l),  is  located  at  axial  station 
12.72,  and  HL(1)  positions  the  hinge  line  2.28  inches  aft  of  the 
root  chord  leading  edge.  The  tip  and  root  chords,  CT(1)  and 
CR( 1 ) ,  are  0.25  and  4.00  inches,  respectively. 

Input  item  6  is  similar  to  item  5,  except  the  geometry  and 
location  of  the  second  fin  set  is  input.  SPAN(2)  =  3.50,  XLE(2) 
=  43.00,  HL (  2 )  =  3.00,  CT(  2)  =  3.50,  and  CR(  2)  =  7.00  define  the 
geometry  of  fin  set  2.  These  values  are  also  specified  in 
inches. 

Input  Item  7  contains  the  NCOND  attitude  conditions  for 
which  results  are  to  be  obtained.  For  this  sample  case,  two  roll 
angles  with  and  without  pitch  control  deflection  are  con- 


sidered.  For  the  first  attitude  condition,  the  body  angle  of 
attack  is  8.0°,  the  roll  angle  is  0.0°,  and  the  deflections  are 
0.0°.  The  second  case  is  the  same  as  one  except  fins  2  and  4  are 
deflected  15°  for  pitch  control.  Case  3  is  the  same  as  case  1 
except  that  the  roll  angle  is  45°,  and  case  4  is  the  same  as  case 
3  except  all  the  fins  are  deflected  15°. 

Input  Item  8  contains  the  run  Mach  numbers.  For  this  case 
there  are  two,  1.30  and  1.75. 

Input  Item  9,  IFLB,  is  a  single  integer  flag  indicating 
whether  the  body  forces  are  to  be  input  or  calculated  by 
MISLE3 .  For  this  case  they  are  input,  IFLB  =  0.  Input  item  10 

'fc 

is  a  set  of  NMACH  NCON  cards  containing  the  axial-force,  normal- 
force,  and  pitching-moment  coefficients  for  the  body  alone.  All 
the  values  are  input  as  zero  for  this  case. 

Output  Description 

A  typical  output  file  from  MISLE3  is  described  in  this  sec¬ 
tion.  In  general,  the  output  quantities  from  MISLE3  are  labeled 
and  each  page  is  headed  with  appropriate  descriptive  informa¬ 
tion.  The  output  from  Sample  Case  1  requires  four  pages  which 
are  described  in  this  section  and  shown  in  Figure  17. 

The  first  output  paae  contains  the  MISSILE  3  header  to  indi¬ 
cate  the  version  number  and  creation  date,  the  summary/title 
cards  for  the  run,  and  an  echo  of  the  input.  The  program  control 
variables,  reference  information,  fin  geometries,  attitude  condi¬ 
tions,  Mach  numbers,  and  body  nose  information  are  all  printed 
for  checkina  purposes.  Along  with  the  echo  of  the  input  fin 
aeometry,  the  calculated  aspect  ratio  and  taper  ratio  are 
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The  results  from  MISLE3  start  on  output  page  2  which  con¬ 
tains  the  output  for  the  first  attitude  condition.  The  run 
number,  Mach  number,  body  angle  of  attack,  and  roll  angle  are 
printed  on  the  second  line.  The  next  block  of  output  contains 
the  component  and  overall  loads  acting  on  the  missile.  The  nose 
loads  include  components  from  the  potential  normal-force  coeffi¬ 
cient  slope,  dCN/da,  and  from  a  crossflow  drag  term.  Fin  section 
1  loads  are  the  sum  of  the  loads  on  the  individual  fins  of  fin 
set  1  along  with  the  body  carryover  loads.  The  afterbody  loads 
are  due  to  the  presence  of  fin  and  body  vortices  over  the  after¬ 
body.  Fin  section  2  loads  are  a  sum  of  the  individual  fin  loads 
and  the  body  carryover  loads.  The  total  loads  on  the  missile  are 
the  sum  of  the  component  loads. 
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The  next  block  of  output  on  page  2  contains  the  individual 
fin  loads  for  both  fin  sets.  The  four  fins  of  fin  set  1  are 
given  first  followed  by  fin  set  2.  The  deflection  angle,  if  non¬ 
zero,  is  printed  followed  by  the  equivalent  angle  of  attack,  the 
normal-force  coefficient,  the  hinge-moment  coefficient,  the  bend¬ 
ing-moment  coefficient,  and  the  rolling  moment  coefficient.  This 
completes  the  output  for  the  first  attitude  condition. 

The  output  for  attitude  conditions  2  and  3  are  contained  on 
pages  3  and  4,  respectively,  and  have  the  same  format  as  output 
page  2.  This  completes  the  output  description. 

RESULTS 

This  section  pLesents  results  obtained  with  the  engineering 
prediction  method  MISSILE  3.  Comparisons  of  the  results  to 
experimental  data  are  made  for  code  verification  and  to  indicate 
necessary  improvements.  This  section  is  divided  into  two  sec¬ 
tions;  Body-Tail  configurations  and  Canard-Body- Ta il  configura- 


Body-Tail  Configurations 


Two  independent  comparisons  with  experimental  data  on  body- 
tail  geometries  are  discussed  in  this  section.  The  first  con¬ 
figuration  is  a  2.5-caliber  ogive  nose  followed  by  10.9  diameter 
cylindrical  body  (Ref.  27).  The  tail  fins,  located  11.2  diame¬ 
ters  aft  of  the  nose  tip,  have  an  aspect  ratio  of  2.81,  a  taper 
ratio  of  0.423,  and  a  body  radius  to  fin  semi-span  ratio  of 
0.25.  This  missile  is  depicted  in  Figure  18. 

Comparisons  of  measured  and  predicted  aerodynamic  character¬ 
istics  at  M  =  1.6  are  shown  in  Figure  19  -  22.  Fiqure  19(a) 
illustrates  the  overall  normal-force  coefficient  and  the  pitch¬ 
ing-moment  coefficient  at  zero  roll  angle  with  no  fin  deflec¬ 
tions.  The  normal  force  is  in  good  agreement  with  experiment  at 
low  and  moderate  angles  of  attack  and  slightly  underpredicted  at 
high  angles  of  attack.  The  pitching  moment  is  also  in  good 
agreement  for  low  and  moderate  angles  of  attack  but  overpredicted 
at  high  angles  of  attack  above  12  degrees.  The  predicted  loads 
on  fin  2  are  compared  with  experiment  in  Figure  19(b).  The  fin 
normal-force  coefficient  comparison  indicates  good  agreement  for 
fin  loads,  but  the  hinge-moment  coefficient  comparison  is  not  as 
good.  Hinge  moments  are  sensitive  to  thickness  distribution,  and 
for  this  case  the  fins  have  a  wedge  leading  edge  then  a  flat 

constant  thickness  section  followed  by  a  wedge  trailing  edge. 

The  data  base  fins  in  MISLE3  are  all  double  wedge  sections.  The 

fin  bending-moment  coefficients  are  in  good  agreement  with  expe¬ 
riment. 

Figures  20(a)  and  20(b)  show  results  for  M  =  1.6  at  zero 

roll  angle  with  fins  2  and  4  deflected  -15  deqrees  for  nose  up 
pitch  control.  The  overall  normal  force  in  Figure  20(a)  and  the 
fin  normal  force  in  Figure  20(b)  are  in  good  agreement  with  mea¬ 
sured  loads;  however,  the  pitching  moment  comparisons  in  Figure 
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20(a)  indicates  an  overprediction  of  complete  configuration 
pitching  moment.  Since  the  fin  normal  force  is  predicted  well, 
the  difference  in  the  pitching  moment  arises  from  an  e^ror  in  the 
body  load  distribution.  The  fin  loads  are  small  due  to  the  nega¬ 
tive  deflection  angle  and  the  positive  angle  of  attack;  there¬ 
fore,  the  error  in  the  body  load  distribution  dominates  the 
pitching  moment.  This  is  not  the  case  when  the  fin  is  highly 
loaded  in  the  undeflected  case  shown  in  Figure  19.  The  error  in 
the  body  load  distribution  arises  from  using  the  crossflow  drag 
formulation  for  the  normal  force  on  the  body.  This  correlation 
predicts  overall  normal  force  well  but  does  not  provide  a  force 
which  varies  along  the  body. 

Figure  21  compares  results  for  M  ra  =  1.6  at  45  degrees  roll 
without  fin  deflections.  Figure  21(a),  (b),  and  (c)  indicate  the 
normal  force  on  the  fins  and  body  are  in  good  agreement  with 
experiment.  The  results  for  the  lee  side  fins  in  Figure  21(b) 
show  that  the  effect  of  the  body  wake  on  the  fins  is  handled  well 
by  the  program  for  this  case.  Figure  21(c)  indicates  an  under¬ 
prediction  of  normal  force  •  n  the  windward  side  fins.  Bending 
moments  are  generally  in  good  agreement  for  this  case. 

Figure  22  compares  results  for  M  =  1.6  at  45  degrees  roll 
with  -15  fin  deflection  on  all  four  fins  for  positive  pitch  con¬ 
trol.  The  results  are  similar  to  the  results  for  pitch  control 
at  zero  roll  angle  in  Figure  20.  The  overall  and  fin  normal 
forces  agree  well  with  experiment,  but  the  pitching  moment  is 
overpredicted.  This  is  attributed  to  a  deficiency  in  the  body 
load  distribution.  It  is  apparent  from  Figures  20(b),  22(b),  and 
22(c)  that  the  effect  of  control  deflection  on  the  fin  bending- 
moment  coefficient  is  in  error  for  this  fin. 


Figures  23  -  26  show  comparisons  of  measured  and  predicted 
results  on  the  same  missile  model  for  a  Mach  number  of  3.7. 
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1.6  results 


Generally  the  results  are  similar  to  the  Mro  = 
described  above.  The  overall  and  fin  normal  forces  are  in  good 
agreement,  and  the  pitching  moment  is  overpredicted  when  fins  are 
deflected  such  that  the  fin  normal  force  is  near  zero  or  when  the 
body  loads  are  large  with  respect  to  the  fin  loads. 

An  overall  view  of  the  body-tail  comparisons  in  Figure  19  - 
26  indicate  that  MISSILE  3  predicts  the  missile  and  fin  loads 
well.  There  is  a  deficiency  in  the  body  load  distribution  model, 
and  the  user  should  input  body  loads  if  a  better  body  model  or  if 
correlations  from  experiment  are  available. 

The  second  body-tail  model  considered  with  MISSILE  3  is 
shown  in  Figure  27  (Ref.  28).  The  body  is  a  1.5-caliber  ogive 
nose  followed  by  a  15-diameter  cylindrical  afterbody.  The  fin 
aspect  ratio  is  2.57,  the  taper  ratio  is  0.4,  and  the  body  radius 
to  fin  semi-span  ratio  is  0.25.  This  investigation  did  not  pro¬ 
vide  any  fin  loads,  but  total  configuration  normal  force  and 
pitching  moment  for  Mach  numbers  of  0.7,  0.9,  1.42,  2.01,  and 

3.08  are  available  for  comparison  purposes. 

The  measured  and  predicted  normal-force  and  pitching-moment 
coefficients  are  compared  in  Figures  28(a)-(e)  for  zero  roll 
angle  and  no  fin  deflections.  The  Mot  =  0.7  and  0.9  results  in 
Figures  28(a)  and  (b)  indicate  that  normal  force  is  predicted 
well  at  low  to  moderate  angles  of  attack,  but  they  are  over¬ 
predicted  at  moderate  to  high  angles.  Since  fin  loads  are  not 
available  for  this  case,  it  is  not  apparent  why  the  predicted 
normal  force  is  too  large.  It  is  possible  that  the  normal  force 
is  sensitive  to  fin  cross  section  geometry  in  the  transonic 
regime.  The  test  fins  have  wedge  leading  edges,  a  flat  constant 
thickness  section,  and  a  wedge  trailing  edge.  The  data  base  fins 
are  double  wedge  sections. 
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Figures  28(c)  through  (e)  illustrate  results  at  supersonic 
Mach  numbers  of  1.42,  2.01,  and  3.08,  respectively.  These 
results  are  similar  to  the  previous  supersonic  body-tail  results 
described  in  Figures  18  -  26.  The  predicted  normal-force  coeffi¬ 
cient  is  in  good  agreement  with  experiment  for  all  cases,  but  the 
pitching  moment  is  overpredicted  at  moderate  to  high  angles  of 
attack.  The  error  in  the  body  load  distribution  is  likely  the 
cause  of  the  pitching  moment  error,  and  the  error  is  even  larger 
for  this  particular  case  since  an  extremely  long  body  is 
involved. 


Canard-Body- Tail  Results 

Comparisons  of  measured  and  predicted  characteristics  on 
three  canard-body-tail  models  are  described  in  this  section.  The 
first  configuration  is  a  canard-controlled  missile,  similar  to 
the  Sidewinder  missile,  presented  in  Reference  25  and  depicted  in 
Figure  29.  The  body  is  a  2.25-caliber  ogive  nose  followed  by  a 
21.9-diameter  cylindrical  body.  The  canard  fins  have  an  aspect 
ratio  of  2.00,  a  taper  ratio  of  0.30,  and  a  body  radius  to  fin 
semi-span  ratio  of  0.226.  The  canard  leading  edge  is  located 
3.73  diameters  aft  of  the  nose  tip.  The  tail  fins  under  consid¬ 
eration  have  an  aspect  ratio  of  1.06,  a  taper  ratio  of  0.58,  and 
a  body  radius  to  fin  semi-span  ratio  of  0.19.  The  tail  fins  are 
located  12.7  diameters  aft  of  the  canard  trailing  edge.  Overall 
loads  on  the  model  at  several  supersonic  Mach  numbers  are  avail¬ 
able,  but  only  results  for  Mot  =  2.5  are  presented  here.  Refer¬ 
ence  25  also  provides  comparisons  to  MISSILE  2A  which  are  also 
included  herein. 

Measured  and  predicted  normal-force  and  pitching-moment 
coefficients  at  M  =  2.5  for  zero  roll  anqle  and  without  control 
deflections  are  shown  in  Figure  30.  MISSILE  3  underpredicts  the 
normal  force  at  high  angles  of  attack;  however,  for  a  greater 


than  23.6  degrees,  the  crossflow  Mach  number  exceeds  unity  and 
the  correction  for  local  dynamic  pressure  and  local  Mach  number 
activate  as  described  in  the  equivalent  angle  of  attack  sec¬ 
tion.  This  correction  tends  to  increase  the  normal  force  to  a 
level  comparable  to  projected  experimental  results.  Since  this 
correction  helps  the  calculation,  an  effort  should  be  made  to 
extend  the  correction  jowl  to  the  critical  crossflow  Mach  number 
for  a  circular  cylinder.  This  corresponds  to  a  crossflow  Mach 
number  of  approximately  0.5  which  occurs  at  an  angle  of  attack  of 
11.54  degrees  in  this  case.  The  correction  also  improves  the 
pitching  moment  results. 

Figure  31  shews  a  comparison  of  control  deflection 
effects.  Canard  fins  2  and  4  are  deflected  +5  degrees  for  posi¬ 
tive  pitch  control.  The  normal-force  coefficient  is  under¬ 
predicted  at  high  angles  of  attack,  but  the  crossflow  correction 
in  the  equivalent  angle  of  attack  formulation  increases  it  when 
the  crossflow  Mach  number  exceeds  unity. 

Figure  32  compares  measured  and  predicted  results  on  a 
rolled  configuration  without  fin  deflection.  As  with  the  pre¬ 

vious  comparisons,  the  normal-force  coefficient  is  under¬ 
predicted.  In  general,  the  MISSILE  2A  results  are  in  better 

agreement  with  experiment  over  the  entire  range  of  angle  of 
attack. 

The  next  model  considered  is  the  Army  Generalized  Missile  of 
Reference  26  shown  in  Figure  33.  This  model  consists  of  a  3- 

caliber  rounded  ogive  nose  and  a  7. ^-diameter  cylindrical  after¬ 
body.  The  canard  fins  have  an  aspect  ratio  of  3.53,  a  taper 

ratio  of  0.0625,  and  a  body-radius  to  fin  semi-span  ratio  of 
0.42.  The  tail  fins  have  an  aspect  ratio  of  1.33,  a  taper  ratio 
of  0.5,  and  a  body  radius  to  fin  semi-span  ratio  of  0.42.  The 
experimental  test  for  this  missile  provides  overall  loads  and 
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detailed  canard  and  tail  fin  loads.  Results  for  a  Mach  number  of 
1.75  are  presented  in  this  report. 

Figure  34(a)  illstrates  the  compariron  of  measured  and  pre¬ 
dicted  overall  normal-force  and  pitching-moment  coefficients  on 
the  model.  The  normal  force  is  slightly  under predicted  at  high 
angles  of  attack,  and  the  pitching  moment  is  overpredicted  in  the 
same  flow  regions.  Figure  34(b)  indicates  the  loads  on  the 
canard  fins  are  underpredicted  for  high  angles  of  attack.  This 
is  due  primarily  to  the  nose  vortex  flow  model  interaction  with 
the  canard  fins.  In  the  body-tail  case  the  nose  and  forebody 
vortex  model  did  not  cause  this  type  of  an  interference  effect 
since  the  vortex  system  is  further  away  from  the  body  and  fins 
because  of  the  long  forebody.  The  close  proximity  of  the  nose 
vortex  and  the  canard  fins  results  in  the  error  in  the  canard 
force.  The  error  is  even  more  apparent  in  the  roll  cases  which 
follow. 

The  tail  fin  loads  predicted  by  MISSILE  3  in  Figure  34(c) 
agree  extremely  well  with  the  experimental  data.  Based  on  the 
above  comparisons,  the  error  in  the  overall  pitching-moment  coef¬ 
ficient  in  Figure  34(c)  must  be  due  primarily  to  an  error  in  the 
body  load  distribution  since  the  tail  fin  loads  are  in  good 
agreement  with  experiment  and  the  canard  fin  loads  are  under¬ 
predicted  . 

Figure  35  illustrates  the  effect  of  canard  control  deflec¬ 
tion  on  the  forces  and  moments  on  the  unrolled  model.  Canard 
fins  2  and  4  are  deflected  to  -fl5  degrees  for  positive  pitch 
control.  These  results  are  similar  to  those  described  for  the 
undeflected  fin  case.  The  normal  force  is  in  good  agreement  with 
experiment  except  at  high  angles  of  attack  where  it  is  slightly 
underpredicted  as  seen  in  Figure  35(a).  The  canard  fin  loads  are 
underpredicted,  Figure  35(b),  as  in  the  undeflected  case  shown  in 
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Figure  34(b).  The  tail  fin  loads  in  Figure  35(c)  are  in  good 
agreement  with  the  measured  results.  In  general,  the  hinge- 
moment  coefficients  are  in  fair  agreement  with  experiment,  and 
the  bending-moment  coefficients  agree  very  well  with  experiment. 

Figure  36  shows  the  effect  on._the  fin  loads  of  rolling  the 
model  45  degrees.  Figure  36(a)  indicates  the  predicted  normal 
force  is  in  good  agreement  with  experiment  except  at  high  angles 
of  attack  where  it  is  underpredicted.  The  lee  side  canard  loads 
shown  in  Figure  36(b)  indicate  the  individual  fir.  normal  forces 
are  dramatically  underpredicted.  This  is  likely  due  to  the  prox¬ 
imity  of  the  nose  vortex  system  to  the  lee  side  canard  fins. 
This  effect  was  seen  in  the  unrolled  case,  but  it  was  not  as 
dramatic.  The  windward  canard  fin  loads  are  seen  in  Figure 
36(c),  and  these  results  indicate  that  the  windward  fin  loads  are 
in  satisfactory  agreement  with  experiment.  Figure  36(d)  illus¬ 
trates  the  lee  side  tail  fin  loads.  The  normal  force  is  pre¬ 
dicted  well  for  low  to  moderate  angles  of  attack,  but  at  higher 
angles  of  attack,  the  normal  force  is  underpredicted  due  to  the 
afterbody  wake.  The  vortex  wake  also  causes  a  large  discrepancy 
in  the  tail  fin  hinge  moment.  The  windward  side  tail  fin  loads 
arc  in  good  agreement  with  experiment  as  seen  in  Figure  36(e). 

The  effects  of  canard  fin  deflection  on  the  loads  for  the 
model  rolled  45  degrees  are  shown  in  Figure  37.  These  results 
are  similar  to  the  undeflected  case  discussed  above. 

The  final  comparison  with  experiment  is  for  the  canard-body- 
tail  model  shown  in  Figure  38  (Fef.  28).  The  body  is  a  1.5- 
caliber  ogive  nose  followed  by  a  15-diameter  cylindrical  body. 
The  canard  fins  have  an  aspect  ratio  of  0.857,  a  taper  ratio  of 
0.4,  and  a  body  radius  to  fin  semi-span  ratio  of  0.25.  The  tail 
fins  have^an  aspect  ratio  of  2.57,  a  taper  ratio  of  0.4,  and  a 
body  radius  to  fin  semi-span  ratio  of  0.25.  Overall  model  loads 
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The  second  objective  was  to  include  the  new  data  base  into  a 
computer  code  to  provide  an  engineering  prediction  method  for  the 
aerodynamic  characteristics  of  typical  missile  configurations 
over  a  wide  range  of  flow  conditions.  Part  of  this  effort 
included  improvements  to  the  equivalent-angle-of-attack  method, 
investigation  of ..  efficient  means  of  manipulating  large  data 
bases,  incorporation  of  correction  methods  for  local  Mach  number 
and  dynamic  pressure  in  regions  of  supersonic  crossflow,  and 
investigation  of  improved  vortex-induced  effects. 

Both  the  above  objectives  were  accomplished  succesfully.  A 
large  body  of  experimental  data  was  collected  from  four  separate 
tunnel  entries,  and  these  data  are  resident  at  the  NASA/Langley 
Research  Center.  The  magnitude  of  the  data  base  is  such  that  it 
was  impossible  to  analyze  and  correlate  all  the  information  for 
this  investigation.  An  engineering  prediction  method  in  the  form 
of  a  code,  MISSILE  3,  was  developed,  and  verification  was  accom¬ 
plished  by  comparison  with  independent  experimental  results.  The 
general  result  from  this  effort  is  that  an  improved  prediction 
method  for  the  normal  force,  hinge  moment,  and  bending  moment  of 
a  wide  range  of  missile  fins  over  a  wide  range  of  flow  conditions 
is  now  available. 
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In  conclusion,  the  extensive  comparisons  of  measured  and 
predicted  aerodynamic  characteristics  presented  in  this  report 
for  verification  purposes  prompt  the  following  observations.  The 
new  code,  MISSILE  3,  is  capable  of  predicting  the  performance 
characteristics  of  typical  missile  configurations  under  an  exten¬ 
sive  range  of  flow  conditions,  and  it  has  application  as  a  pre¬ 
liminary  design  method.  Generally,  predicted  individual  fin 
loads  are  in  good  agreement  with  experiment,  and  loading  distri¬ 
butions  are  adeauately  predicted  for  preliminary  design  pur¬ 
poses.  Fortunately,  this  was  the  prime  objective  for  the  new 
MISSILE  3  code.  The  simplified  body  load  prediction  method  is 


for  Mach  numbers  0.7,  0.9,  1.42,  2.01,  and  3.08  are  available. 
As  with  the  body-tail  results  from  this  same  reference,  no  indi¬ 
vidual  fin  loads  were  measured. 

Figures  39(a)  through  (e)  illstrate  comparisons  of  measured 
and  predicted  overall  normal-force  and  pitching-moment  coeffi¬ 
cients  for  the  unrolled  missile  for  the  Mach  number  range.  The 
normal  force  is  usually  slightly  underpredicted  with  respect  to 
the  experimental  results.  The  crossflow  Mach  number  correction 
in  the  Mm  =  3.08  condition  in  Figure  39(e)  raises  the  value  of 
the  normal  force  and  improves  agreement  with  experiment.  The 
pitching  moment  is  not  predicted  well  for  transonic  Mach  numbers, 
but  it  exhibits  trends  similar  to  the  other  canard-body-tail 
results  for  supersonic  Mach  numbers. 

Figures  40(a)  through  (e)  contain  the  comparisons  for 
normal-force  and  pitching-moment  coefficients  for  the  above  model 
at  a  roll  angle  of  45  degrees.  The  predicted  normal  force  is  in 
good  agreement  with  experiment  in  all  cases.  The  pitching  moment 
is  predicted  well  in  the  low  and  moderate  angles  of  attack  range 
except  for  the  M  =  0.9  case. 

CONCLUSIONS 

The  original  intended  purpose  of  the  multiphase  program 
described  in  this  final  report  was  twofold.  The  first  objective 
was  to  obtain  an  extensive  missile-fin  data  base  which  would 
permit  development  of  a  broadly  applicable  engineering  predictive 
method  for  calculating  the  aerodynamic  characteristics  of  anard 
cruciform  missiles.  Specifically,  a  systematic  data  base  which 
extended  (1)  the  geometric  fin  characteristics  of  aspect  ratio 
and  taper  ratio  and  (2)  the  flow  characteristics  of  Mach  number, 
angle  of  attack,  deflection  angle,  and  roll  angle  beyond  the 
current  limits  of  existing  experimental  data  was  to  be  developed. 
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satisfactory  for  total  normal  force,  but  the  distribution  of 
normal  force  on  the  body,  dominated  by  vortex  shedding  character¬ 
istics  at  high  angles  of  attack,  is  marginal.  For  this  reason, 
the  provision  for  the  input  of  improved  body  characteristics  is 
included  in  the  code. 

The  performance  of  MISSILE  3  is  very  similar  to  that  of  the 
previous  MISSILE  2A  code,  but  the  new  code  is  applicable  over  a 
much  wider  range  of  geometric  and  flow  parameters.  The  code  is 
simple  and  economic  to  use,  and  the  enclosed  user's  manual  with 
sample  input  and  output  should  expedite  the  learning  process 
associated  with  a  new  code.  However,  it  should  be  noted  by  pros¬ 
pective  users  that  MISSILE  3  is  a  newly  developed  code,  not  a 
revised  version  of  MISSILE  2A.  There  are  likely  to  be  subtle 
errors  that  were  not  uncovered  during  the  extensive  development 
phase  of  the  code;  therefore,  the  user  is  cautioned  to  carefully 
test  predicted  results  against  common  sense  and  experience  before 
accepting  the  results.  The  authors  encourage  both  positive  and 
negative  feedback  from  users  of  the  code. 

As  noted  in  various  parts  of  this  report,  there  are  certain 
limitations  of  the  code  and  suggested  improvements  that  were 
beyond  the  scope  of  the  present  effort.  The  next  section  briefly 
outlines  some  of  these  areas. 

RECOMMENDATIONS 

The  first  and  most  important  recommendation  is  for  addi¬ 
tional  testing  and  verification  of  the  new  MISSILE  3  code  to 
better  define  the  limits  of  its  capability.  Extensive  compar¬ 
isons  with  independent  experimental  results  for  a  wide  range  of 
geometries  and  flow  conditions  will  help  identify  the  practical 
range  of  application  of  the  code  and  uncover  heretofore  unde¬ 
tected  errors. 
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Several  known  areas  of  needed  improvement  include  the  following: 

1.  Modify  the  body  force  distribution  calculation  to 
include  enhanced  vortex  shedding  characteristics  either 
through  empirical  means  or  the  inclusion  of  now- 
available  discrete  vortex -methods  (Refs.  13  and  14). 

2.  Improve  the  fin  hinge  moment  pre  .ction  by  including 
additional  empirical  information  J/or  correlations  for 
the  effects  of  different  airfoil  ;  ictions. 

3.  Extend  the  correction  for  local  Mach  number  and  dynamic 
pressure  effects  to  include  crossflow  Mach  numbers 
between  0.5  and  1.0. 

4.  Extend  the  correlations  for  fin  control  effects  through 
the  incorporation  of  available  control  data  not  con¬ 
sidered  in  the  current  effort. 

5.  Extend  the  method  to  include  effects  of  changing  missile 
afterbody  radius. 

6.  Extend  the  method  to  include  effects  of  noncircular 
cross  section  missiles. 
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LIST  OF  SYMBOLS 


a  body  radius 

AR  aspect  ratio  of  wing-alone  formed  by  joining  two  fins  at 

their  root  chords 

Cft  axial-force  coefficient?  axial  force/qraSR 

CBM  fin  bending-moment  coefficient?  bending  moment/q^Sj^lp 

CHM  fin  hinge-moment  coefficient?  hinge  moment/  qOTSRlR 

rolling-moment  coefficient?  rolling  mcment/q^Sj^lj^ 

CMx  rolling-moment 

CMy  pitching-moment 

CMz  yawing-moment 

Cm  pitching-moment  coefficient?  pitching  moment/q^Sj^ 

CN  normal-force  coefficient?  normal  force/qroSR 

CNW  wing-alone  normal-force  coefficient?  normal  force/g^Sp 

CR  fin  root  chord 

Cx  force  acting  along  the  x-axis,  = 

Cy  force  acting  along  the  y-axis,  side  force 

Cz  force  acting  along  the  z-axis,  =  CN 

D  maximum  body  diameter 

K„  Beskin  upwash  factor 

1R  reference  length 

Mi  average  local  Mach  number 

M  freestream  Mach  number 

oo 

q^  average  local  dynamic  pressure 

am  freestream  dynamic  pressure 

s  exposed  fin  semispan 

sm  semispan  of  fin-body  combination 

SR  reference  area 

v,w  velocity  components  in  the  real  plane 

W  complex  potential 

x  axis  along  the  missile  centerline  positive  aft 

chordwise  location  of  fin  hinge  line  measured  from  lead¬ 
ing  edge  of  the  fin  root  chord 
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XLE 

"x 
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c 
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eq 

6 

aadj 

AOPP 

aself 

X 

XADJ 

XOPP 

4> 


axial  location  of  the  leading  edge  of  the  root  chord 
measured  from  the  nose  tip 

axial  location  of  the  fin  normal  force  center-of-pres- 
sure  measured  from  the  leading  edge  of  the  fin  root 
chord 

axis  to  the  right  looking  forward 

spanwise  location  of  the  fin  normal  force  center-of- 
pressure  measured  outboard  from  the  fin  root  chord 
axis  in  the  wind  plane  positive  up 

included  angle  of  attack;  angle:  between  the  body  axis 
and  the  freestream  velocity  vector. 

equivalent  angle  of  attack;  that  angle  of  attack  of  che 
wing-alone  for  which  its  normal  force  is  twice  that  of 
the  fin 

fin  deflection  angle 

slender-body-theory  factor  for  carryover  of  loading  due 
to  deflection  to  an  adjacent  fin 

slender-body-theory  factor  for  carryover  of  loading  due 
to  deflection  to  an  opposite  fin 

slender-body-theory  factor  for  effectiveness  of  fin 
de  flection 

taper  ratio;  ratio  of  fin  tip  chord  to  root  chord 
fraction  of  the  fin  planform  area  affected  by  the  load¬ 
ing  on  the  adjacent  fin 

fraction  of  the  fin  planform  area  affected  by  the  load¬ 
ing  on  the  opposite  fin 
roll  angle 


Subscripts : 


i  fin  position  numbered  clockwise  looking  forward;  fin  1 

of  the  first  fin  section  is  always  in  the  first  quadrant 
i,0  data  base  value  for  the  i-th  fin 
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1.1  scaled  i~th  fin  value  for  dimension  of  interest  with  no 
fin  deflection 

1.2  scaled  i-th  fin  value  for  dimension  of  interest  with  fin 
deflection  included 

Superscripts : 

A  indicates  a  quantity  which  lacks  the  effects  of  vorti- 

city  in  the  field 
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igure  2.-  Banked  canard-cruciform  missile  at  angle  of  attack 
showing  typical  vortex  field. 


(a)  Total  loads  and  fin  numbering  system. 
Figure  3.  "Symbols  and  conventions. 
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(b)  Sign  convention  for  canard  deflection  angles 


Figure  3.—  Continued. 
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Figure  4.-  Concluded. 
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Figure  5.- 


Stability  normal-force  coefficient  for 
fin  52;  AR  =  2,  A  =  0.5. 


Separation  angle,  0^ ,  degrees 


Landrum,  LRC 
(ref.  18) 
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Figure  11.-  Subroutine  calling  sequence  for 
program  MISLE3. 
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Figure  12.-  Continued. 
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Figure  13.-  Common  block  cross  reference  list  for  program 
MISLE3 . 
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Figure  18.-  Sample  case  2  input. 
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Figure  18.-  Geometry  for  the  body-ta 


Figure  19 .  -  Concluded. 
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□  Experiment,  Ref.  2? 
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Figure  21.- 


(a)  Overall  loads 

Comparison  of  predicted  and  measured 
aerodynamic  characteristics  for  the 
body-tail  model  in  Reference  27 , 

M  =  1.6,  <t>  =  45°,  6.  =  0°. 
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Figure  24.- 

Comparison  of  predicted  and  measured 
aerodynamic  characteristics  for  the 
body-tail  in  Reference  27,  M  =3.7, 
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□  Experiment,  Ref.  27 
—  MISSILE  3 
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Figure  25.-  Concluded. 
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MISSILE  3 
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Figure  30.-  Comparison  of  predicted  and  measured 
aerodynamic  characteristics  for  the 
canard-body-tail  model  in  Reference  25, 
M  =  2.5,  <p  =  0°,  6.  =  0°. 
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APPENDIX  A 


Local  Mach  Number  and  Dynamic  Pressure  Correction 

A  body  at  angle  of  attack  in  which  the  crossflow  Mach  number/ 
Moo  sin  aC/  exceeds  unity  experiences  regions  of  flow  where  local 
conditions  are  very  different  from  freestream  conditions..  For 
example,  the  presence  of  a  bow  shock  wave  can  significantly  alter 
local  Mach  numbers  and  dynamic  pressures;  therefore,  a  method 
incorporating  the  equivalent  angle  of  attack  approach  which 

relies  on  wing  alone  aerodynamic  characteristics  in  a  uniform 
freestream  requires  a  correction  of  the  referenced  flow 

conditions  to  account  for  the  local  flow  conditions  seen  by  the 

fins  on  the  body.  More  detailed  descriptions  of  this  phenomena 
are  presented  in  References  A1  and  A2. 

A  correction  for  local  dynamic  pressure  and  Mach  number 
effects  is  implemented  using  predicted  velocity  components, 
densities,  ar  pressures  from  an  Euler  code.  The  velocity 

components  u,  v,  and  w,  the  density,  p,  and  the  pressure,  p,  are 
determined  around  a  body  using  program  SWINT.  The  body  consists 
of  a  3-caliber  ogive  nose  followed  by  a  long  cylindrical  section, 
and  the  flow  field  properties  are  investigated  10  diameters  afc 
of  the  nose  section.  The  following  dimensionless  variables  are 
defined 

dimensions  x  =  x/1  ,  r  =  r/1  ,  <p  =  f 

axial  velocity  w  =  77-— - 

U  cosa 
00 

(Al)  Hemsch,  M.  J.  and  Nielsen  J.  N. :  Triservice  Program  for 
Extending  Missile  Aerodynamics  Data  Base  and  Prediction 
Program  Using  Rational  Modeling.  Interim  Report  for  Period 
Juno  1?.,  1982  to  June  15,  1983.  NEAR  TR  305,  Aug.,  1983. 


(Al) 

(A2) 


lateral  velocities 

~  u 

U  =  tt  • 

U  sma 

04 

(A3) 

~  V 

U  sina 

00 

(A4) 

density 

p  =  p/p«, 

(A5) 

pressure 

~  .  p 

P  “  „2  .  2 

Poo13*32-*  a 

(A6) 

Note  that  a,  V,  Ht,  p,  and  p  are  equal  at  corresponding  points  in 
the  flow  field  (X,f,$)  if  the  parameters  ki  =  tfcota  and  k2  = 
Moosina  have  the  same  valves  at  the  points.  For  a  long  (infinite) 
swept  cylinder,  k^  -*■  0,  X  <*>,  and  all  of  the  dimensionless 
variables  become  functions  of  Moosina  only. 

Considering  the  pressure  equation,  since  U&  =  M^7poo//?oo  for 
isentropic  flow.  Equation  (A6)  can  be  written  as 


*  -  w2  :  2  <A7> 

or 

^  p  w2  *2  /AON 

P  =  ~  =  7P  Moos:>-n  a  (A8) 

^00 

Since  p  =  p(Moosina)  for  ki  0,  X  <»,  Equation  (A8)  has  the  form 


P*  =  =  f(Mwsin  a)  (A9) 

MX) 


(A2)  Hemsch,  M.  J. ,  and  Mullen,  J. ,  Jr.:  Analytical  Extension 
of  the  MISSILE  1  and  MISSILE  2  Computer  Programs,  March 
1982,  NEAR  TR  272. 
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LiV 
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The  normalized  dynamic  pressure  is 


% 


='M2+  £)  ♦  JfcM 


(A10) 


Substituting  Equations  (A2)  -  (A5)  into  (A10)  yields 


q  ~  T  .  2  ,~2  ~2.  2  ~2.  I 

--a—  =  p  sm  a(u  +  v  )  +  cos  a(w  ) 

%>  L 


(All) 


The  average  local  dynamic  pressure  at  a  given  <p  is  obtain  by 
integrating  Equation  (All)  over  the  exposed  span  and  dividing  by 
the  exposed  span  length;  therefore. 


.  2 

q  _  sm  a 
%>  ~ 


s 

,m 


s  -  a 
m 


2 

~/~2,  ~2nj,  ,  cos  a 

p( u  +  v  )dt  + 


s 


s  -  a 
m 


p (w2) dt 


(A12) 


The  normalized  Mach  number  can  be  expressed  as 


,,  .2  2  2  2  7p  Ip 

M  )  u  +  v  +  w  .  /poo'  ^oo 


ot)  ■ 


IV /p 


U2 

00 


p/p 


00 


p/p, 


J(t) +  (fc) +  ] 


(A13) 


Substituting  Equations  (A2)  -  (A5)  into  (A13)  yields 


fM  )2  ~p\  .  2  ,~2  ~2.  2  /~2v  I 

liiTJ  =  p*isin  a(u  +  v  )  +  cos  a(w  )J 


(A14) 


00 


The  average  local  squared  Mach  number  is 


pq  = 


2  .  2 
sm  a 


s 


e  -  a 
m 


~  o  o  2 

/9  ,-2 .  ~2S  ,,  .  cos  a 

(u  +  v  )dt  + - 

Jr 


s  -  a 
m 


s 

£*(w2)dt 

a 


(A15) 


Note  that  Equation  (A15)  gives  the  average  M2  not  the  square  of 
M. 


-A3 


jr  v*  Vu  ^  -•  -  +v  A.  **  *•  -«Sr  <Ar  A  AV.  *V  Aj  ^  A-  A  «A  A  -\  A.'  A  A;  "A-  < :  A>  ’  A  A  -avv  A  A/  *Ar 


The  following  coefficients  are  defined  for  use  of  the  SWINT- 
generated  data  base  at  x/D  =  10. 


ACr^M^sina)  5  />(u2+  v2) 

(A16) 

BCr^M^sina)  =  /jw2 

(A17) 

c(r>0>Mcosina)  H  A/p* 

(A18) 

=  B/p* 

(A19) 

The  average  values  integrated  over  the  fin  span  are 


A(a/sm,^,MMsina)  =  — 


ItJa  (i',^,M<x)sina)dr  (A20) 


m  a 
s 


B(a/sm,^,MMsinffi)  =  (r  >  1> ,  M^sina)  dr  (A21) 

m  a 

Sm 

C(a/sm,^,Moosina)  5  (r , i , M^sina) dr  (A22) 


m  a 

sm 

D(a/sm,  ^M^sina)  =  (r ,  <t> ,  MoasinQ:)  dr  (A23) 

m  a. 

Using  Equations  (A20)  -  (A23),  the  local  dynamic  pressure  and 

Mach  number  become 


q  r  .  2  ^  2 

—  asiu  a  +  ccos  a 

% 


rw  2 


fM  ]  ~  fM  I  p;  .  2  p:  2 

]  =  =  ^sin  a  +  Dcos  a 


(A24) 


(A25) 


Tables  of  A,  B,  C,  and  D  were  created  from  SWINT  results  at 
x/d  =  10  and  the  a/sm  in  the  data  base  for  various  crossflow  Mach 
numbers  greater  than  1.0.  These  tables  are  included  in  MISSILE 
3,  and  the  average  local  dynamic  pressure  and  Mach  number  are 
calculated  from  Equations  (A24)  and  (A25),  respectively. 
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